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ABSTRACT 

We study large-scale outflows in a sample of 96 star-forming galaxies at 1 < z < 2, using near-UV 
spectroscopy of Fe II and Mg II absorption and emission. The average blueshift of the Fe II interstellar 
absorption lines with respect to the systemic velocity is —85 ± 10 km s _1 at z ~ 1.5, with standard 
deviation 87 km s _1 ; this is a decrease of a factor of two from the average blueshift measured for far-UV 
interstellar absorption lines in similarly selected galaxies at z ~ 2. The profiles of the Mg II AA2796, 
2803 lines show much more variety than the Fe II profiles, which are always seen in absorption; Mg II 
ranges from strong emission to pure absorption, with emission more common in galaxies with blue UV 
slopes and at lower stellar masses. Outflow velocities, as traced by the centroids and maximum extent 
of the absorption lines, increase with increasing stellar mass with 2-3er significance, in agreement with 
previous results. We study fine structure emission from Fe II*, finding several lines of evidence in 
support of the model in which this emission is generated by the re-emission of continuum photons 
absorbed in the Fe II resonance transitions in outflowing gas. In contrast, photoionization models 
indicate that Mg II emission arises from the resonant scattering of photons produced in H II regions, 
accounting for the differing profiles of the Mg II and Fe II lines. A comparison of the strengths of the 
Fe II absorption and Fe II* emission lines indicates that massive galaxies have more extended outflows 
and/or greater extinction, while two-dimensional composite spectra indicate that emission from the 
outflow is stronger at a radius of ~10 kpc in high mass galaxies than in low mass galaxies. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift 



1. INTRODUCTION 

The evolution of galaxies cannot be understood with- 
out the simultaneous consideration of the surrounding 
gas. Galaxies drive powerful, galactic-scale outflows, 
and must also accrete gas from the intergalactic medium 
(IGM) . Thus the process of galaxy evolution reflects the 
cycling of baryons in and out of galaxies. 

The importance of galactic outflows is supported 
by a wide range of models and observations. An 
incomplete list of phenomena that may be ex- 
plained by galactic outflows includes the shape of 
the mass-metallicity re lation and th e metal enrich- 
ment of the IGM (e.g. iGarnettl l2002t ITremonti et all 
2001 lOppenheimer fe Pavel 120061: lErb et all 12 006a; 
Dalcantonl l2007t iFinlator fe Pavel 120071: iBrooks et ail 



2007t iMannucri et al.l 120091 : iPeeples fe Shankarl 120111 : 
Menard et alJ 1201 If ) . the quen ching of star f orma- 
tion in massi ve galaxies Ce.g. ITremonti et al.l 120071 : 
iHopkins etal] 120081 : iGabor et all 120111 ). and the mis- 
match between the galaxy luminosity fun ction and the 
mass function of dark matter halos (e.g. iBensqn et alJ 
l200l lKeres"et~aT1[200l lOppenheimer et alJl2010j ). Pow- 
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erful galactic outflows must also modulate the in- 
flow of gas into galaxies, constraining the assembly 
of the ba ryonic component and regulating star forma- 
tion (e.g. iPave et al.ll2011t iFaucher-Giguere et alj|2011t 
Ivan de Voort et alJl2012TT 

In the broadest terms, galactic outflows are powered 
by energy and momentum injected into the interstel- 
lar medium (ISM) of galaxies by star formation and/or 
active galactic nuclei (AGN). Hot (~ 10 7 K), metal- 
enriched supernova ejecta may be expelled from galax- 
ies, entraining significant amounts of the cool (~ 10 K) 
ISM in the process. Even before massive stars explode 
as supernovae, radiatio n pressure may acce lerate the cool 
gas to large velocities (|Murrav et alj|2lnTl ) . The detailed 
physics of the outflows remains controversial, however, 
and observational constraints on the velocity of outflows 
as a function of radius, the mass outflow rate, the spa- 
tial extent of outflows, and the scaling of these quantities 
with galaxy mass and star formation rate are needed. 

Galactic winds are most commonly traced via the 
blueshifts of interstellar absorption lines from cool, out- 
flowing gas, seen against the backlight of the stellar 
continuum. In the local universe, outflows are found 
in st arburst galaxies, from dwarf starbursts to ULIRGs 
(e.g. IHeckman et all 119901: iStrickland fe Stevens! 120001: 
i Veilleux et alj|2005t iHenrv et al.ll20Q7t iMartin fe Bouchel 
120091 : ISoto et al.ll2012t ). Both observations and models 
suggest a critical star formation rate density in order 
to drive outflows (IHeckman! 120021 : iMurrav et all 120111 : 
iKornei et alJ 120121 ). and outflow velocity is observed 
to increase with galaxy mass a nd star formation rate 
(|Martinll200l iRupke et al.ll2005T) . 

At higher redshifts, 1.5 < z < 4, outflows are ubiq- 



Fell UVI 



ERB, QUIDER, HENRY AND MARTIN 

Fell UV2 



Fell UV3 



. Z 6 D° 3 /1 



3dHp 



a 6 Di/2 

5/2 
7/2 
9/2 

3d«4s 



9/2 



11/2 
3d 6 4p 



a 6 Di/2 
in 

S/2 
7/2 
9/2 

3dMs 



.z*PV 

4/2 
7/2 

3d 6 4p 



a s D|fl 

3/2 
S/2 

7/2 
9/2 



3d 6 4s 



Fig. 1. — The energy level diagrams for selected transitions of the Fe II atom (based on Figure 7 of Hartigan et al. 1999 and Figure 1 of 
Prochaska ct al. 2011). Each transition is labeled with its rest wavelength (A) and Einstein A-coefncient (s~ 1 ). Upward arrows indicate 
resonance absorption (absorption from the ground state). The fine structure splitting of the energy levels allows for fine structure emission. 
The downward arrows show the Fe II* (fine structure) emission that is energetically tied to the resonance absorption transitions. Only 
Fe II* transitions linked to the resonant absorption transitions are shown. 



uitous, with most galaxies observed show ing absorption 
lines blueshifted by ~ 100-300 km s" 1 (e.g.lShaplev et al 



2003; Weiner ct ai. 2009; Stcidcl et al. 20101: 1 Jones eTal 



20111) . Outflows in high redshift galaxies are also traced 



by Lya, which is seen in redshifted emission as well as 
blueshifted absorption; the emission is interpreted as res- 
onant scattering of Lya photons from receding gas on the 
far side of the galaxy. Because the high redshift studies 
have been restricted to a smaller range in mass and star 
formation rate than the local samples (which extend to 
low mass dwarfs) , trends between outflow properties and 
galaxy properties have been more difficult to determine; 
however, recent results at z ~ 1.4 suggest simila r scalings 
with galaxy mass and SFR (IWeiner et al.ll200l . 

A major difficulty in the study of outflows is the mea- 
surement of their spatial extent, since absorption line 
studies of galaxies offer no information on the location 
of the absorbing gas. Spatial information can be ob- 
tained by measuring absorption in the spectra of back- 
ground objects (QSOs or galaxies) with li nes of sight 
passing near galaxies. iSteidel et al.l (|2010D used com- 
posite spectra of pairs of foreground and background 
galaxies, stacked as a function of impact parameter, to 
show that significant metal line absorption is seen around 
galaxies to distances of ~ 60 kpc, while Lya absorption 
extends to ~ 250 kpc. Searches for star-forming coun- 
terparts to Mg II absorption systems also suggest that 
galax i es drive outflows to large distances (|Bouche et al.l 
120071 IMenard etahl 120111 ). while the observation that 
strong Mg II absorption systems are preferentially found 
along the minor axis of disk galaxies is further evi- 
dence that these systems are associated wi th outflows 
(|Kacprzak et al.l I2012t IBouche et all 1201 2aD . On the 
other hand, n ot all absorption sys tems have detectable 
counterparts (jBouche et al.l l2012bl) and the absorption 
has also been inte rpreted as being due to infalling gas 
(jChen et al.l 120101) . a scenario potentially supported by 
the fact that an increase in Mg II absorp tion is also seen 
along the major axes of disk galaxies feacprza k et al.l 
l20ia IBouche et alJl20r2al) . 

The extent of outflowing gas can also be measured if 
galactic winds can be observed in emission. Extended 



Ha and X-ray emis sion is observed in local starbursts 
(jLehnert et ai.lll999D . but such emission becomes unde- 
tectable at higher redshifts. At z = 0.69, iRubin et al.l 
(|2011|) detect Mg II emis sion at a radius of ~ 7 kpc from 
a starburst galaxy, while ISteidel et al.l (|2011[ ) see diffuse 
Lya halos extending to ~ 80 kpc around star-forming 
galaxies at z ~ 2-3, using stacked narrow-band images. 
Both the Mg II and the Lya photons are interpreted as 
arising from resonant scattering in outflowing gas. 

In this study we focus on Mg II and Fe II absorp- 
tion and emission in galaxies at 1 < z < 2. The 
Mg II AA 2796, 2803 doublet is a particularly useful 
diagnostic because it can be observed from the ground 
over a wide range in redshift, from z ~ 0.2 to z ~ 2, 
thus in principle offering direct comparisons of outflow 
properties over ~ 8 Gyr of cosmic time. In prac- 
tice, samples across this redshift range have so far dif- 
fered significantly in size, selection technique, and type 
of ga l axy observed dTremonti et aLl I2007t IWeiner et all 
20091: IMartin fe Bouchel 120091 IRubin et al.l 120101 | 2011 [ 



Coil et all 120111: iKornei et all l2012t IMartin et all 12012; 
Diamond- Stanic et al. 1 120121) . 

The study of Mg II and Fe II absorption and emission 
has also benefited from r ecent radiative tr a nsfer model- 
ing of these transitions. iProchaska et all (|2011l ) use a 
variety of different galactic outflow models to predict the 
strengths of these lines, emphasizing the importance of 
photon scattering and the effect of re-emitted photons 
on the line profiles. Photons absorbed in resonance lines 
such as Mg II are re-emitted near the systemic velocity, 
creating underlying emission which fills in the absorption 
and may affect inferences about the velocity and cover- 
ing fraction of outflowing gas. Because most of the Fe II 
resonance lines are coupled to fine structure transitions, 
as shown in the energy level diagram in Figure [I] the 
photons absorbed in the Fe II transitions may instead be 
re-emitted in the Fe II* fine structure lines. This paper 
focuses on the effect of photon scattering on measure- 
ments of absorption line velocity centroids and equivalent 
widths, and on Fe II* and Mg II emission and absorption 
as tracers of outflowing gas. 

The paper is organized as follows. We describe our 
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sample and the acquisition and analysis of the data in 
Section [2[ the determination of systemic redshifts in Sec- 
tion [3[ and the determination of stellar population prop- 
erties in Section |4l Results from individual spectra are 
reported in Section [5j in Section 15.11 we discuss Mg II 
emission, and Section 15.21 compares the velocities and 
equivalent widths of the Fe II and Mg II transitions. In 
Section|6]we construct composite spectra based on stellar 
population properties (Section I6.ip and on the presence 
or absence of Mg II emission. We discuss the velocity 
centroids and equivalent widths of the transitions in Sec- 
tions 16.2.11 and 16.2.21 and the maximum observed out- 
flow velocities in Section [6.2.3l The remainder of the pa- 
per, Section 16731 is devoted to the examination of Fe II* 
and Mg II emission. In Section [7] we summarize our re- 
sults. We use a cosmology with Hq = 70 km s _1 Mpc -1 , 
Q m = 0.3, and = 0.7. 

2. OBSERVATIONS AND DATA REDUCTION 

The galaxies observed in this study were d rawn from 
the survey described by iSteidel et al.l ([2004) . Targets 
were photometrically pre-selected using the rest-frame 
UV "BM" and "BX" color critieria; the BM criteria tar- 
get galaxies in the redshift range 1.4 < z < 2, and the 
BX criteria select galaxies primarily at 2 < z < 2.5. 
We observed galaxies in the Q1623 (a — 16:25:45, 6 = 
+26:47:23), Q1700 (a = 17:01:01, 8 = +64:11:58) and 
Q2343 (a = 23:46:05, 8 = +12:49:12) survey fields. 

Observations were conducted with the DEep Imag- 
ing M ulti-Ob ject Spectrograph (DEIMOS; iFaber et all 
120031 ) on the Keck II telescope in August 2008 and 
August 2009. DEIMOS is a medium resolution opti- 
cal spectrometer with spectral coverage from ~ 4100 — 
11,000 A; this range includes rest- frame near-UV features 
from Fe II A 2260 to [O II] AA 3727, 3729 in galaxies with 
1 < z < 2. With the use of a slitmask, the 16.7' x 5.0' 
field of view allows the simultaneous observation of up 
to ~ 100 galaxies. 

We observed one mask in each of the three fields, 
covering a total of 147 galaxies. Galaxies were se- 
lected for the masks based on the following priorities: 1) 
k nown redshift (from the previous far-UV spectroscopy 
of ISteidel etaLll2004D in the range 1.4 < z < 2; 2) the 
BM (as opposed to BX) color criteria, which target the 
redshift range for which features of interest fall within the 
spectral co verage; and 3) brighter galaxies with 1Z < 24. 
We used the 600 1/mm grating, which gives a dispersion 
of 0.65 A/pixel. Exposure times ranged from 7.5 to 11 
hrs, with average seeing FWHM ~ 0.6". Details of the 
observations are given in Table Q] 

Most of the data reduction was done with the 
DEEP2 data processing pipeline ([Newman "eTaTl 120121 : 
iCooper et al.l [2012)0 which produces one-dimensional 
air wavelength-calibrated galaxy, inverse variance and 
sky spectra. We then removed the instrumental signature 
from the data by putting the galaxy and error spectra on 
a relative flux scale with the IRAF software package, us- 
ing standard stars observed in the same manner as the 
galaxies. 

We also removed the strong telluric absorption bands 
due to molecular oxygen at ~ 6800 A and ~ 7500 A. This 

5 http: //www2 .keck.hawaii . edu/inst/deimos/pipeline . html 



correction is usually done by scaling the strength of the 
telluric absorption seen in a standard star spectrum to 
the absorption strength seen in an object's spectrum. 
Because of the long integration times used in this study, 
however, the galaxies were observed over a range of air- 
masses, making a single standard star observation taken 
at a single air mass a poor match to the telluric fea- 
ture strength in the individual galaxy spectra. A better 
measurement of the shape of the telluric features can be 
obtained from an average observed-frame galaxy spec- 
trum for each field (because the galaxies lie at different 
redshifts, combining in the observed frame averages over 
the galaxy properties but gives a good measurement of 
the sky features). We used these average spectra to cre- 
ate telluric feature templates for each field. Each average 
galaxy spectrum was normalized and smoothed, and then 
portions of the spectrum which were outside of the tel- 
luric bands were forced to have a normalized flux value 
of one, leaving the telluric bands as the only deviations 
from continuum. The resulting spectrum served as the 
template telluric spectrum for each field. This template 
was then divided into each of the galaxy and error spec- 
tra. 

Finally, the individual spectra were converted to a vac- 
uum wavelength scale and rebinned by a factor of two to 
increase S/N; the 600 1/mm grating is DEIMOS' lowest 
resolution mode and so results in oversampled data. The 
final pixel scale is 1.25 A per pixel in the observed frame, 
and the spectral resolution (measured from the widths of 
the night sky lines at ~ 7500 A) is ~ 3.7 A (FWHM) or 
~ 150 km s" 1 in the Q1623 field, which was observed 
with 1" slits, and ~ 4.5 A or ~ 175 km s _1 in the Q1700 
and Q2343 fields, which were observed with 1.2" slits. 

We measured redshifts z > 1 for 96 of the 147 galaxies 
observed, either from [O II] AA3727, 3729 emission or Fe II 
and Mg II absorption as discussed in Section [3] below 
(the remaining 51 galaxies were mostly among the fainter 
objects on the masks, and most had spectra too noisy to 
enable redshift determination; 4 galaxies had z < 1). 
The analysis presented here is based on this sample of 96 
galaxies. 

After measuring redshifts, the 96 sp ectra were n orma l- 
ized using the procedure outlined by iRix et afl ([2004) . 
Normalization was done only in the rest-frame wave- 
length range 2200-3000 A; this range contains the Fe II 
and Mg II features of interest, and increasing sky noise at 
redder wavelengths makes normalization more difficult. 
We shifted the spectra into the rest frame, identified re- 
gions of the spectrum which contained flux solely from 
the continuum, and fit a spline curve to both the mean 
wavelength and mean flux values within these defined 
continuum windows. The same windows were used for 
each galaxy, with adjustment on an individual basis to 
avoid the gap between detectors or noise spikes. Each 
galaxy spectrum was then divided by the spline fit to 
its continuum. The uncertainty in the continuum fits is 
~ ±5 % for most galaxies. These continuum normalized 
spectra were used for the remainder of the analysis. 

We also make use of multi- wavelength imaging in order 
to model the spectral energy distributions of the galax- 
ies. The opti cal images, in the U „,, G and 1Z filters, are 
described by ISteidel etahl (120041) . the near-IR (J and 
K s ) images bv lErb et al.l (120061)1) . and the mid-IR Spitzer 
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TABLE 1 

Observations 



Field Date Exposure Time ^central 

Slit Width Blocking Filter 
(UT) (s) (A) (arcsec) 

Q1623 02-04 Aug 2008 34,200 7500 1.0 OG550 

Q1700 22-24 Aug 2009 27,000 7200 1.2 GG495 

Q2343 02 Aug 2008; 21-24 Aug 2009 39,600 7200 1.2 OG550 



IRAC data by IShaplev et ail (|2005D and IReddv et all 
(f20f^l2012h . 

3. DETERMINATION OF SYSTEMIC REDSHIFT 

The systemic redshift of a galaxy (i.e., the redshift of 
its stars) serves as the velocity zero-point for studying 
the kinematics of galactic and circumgalactic gas, and 
is used to align galaxies in the rest-frame for the con- 
struction of composite spectra, as described in Section 
O A reliable determination of the systemic redshift is 
therefore required for this study. 

The preferred method of redshift determination, mea- 
surement of the wavelengths of strong nebular emission 
lines from H II regions, becomes increasingly difficult for 
galaxies with z > 1.5, as these lines redshift into the near- 
IR. The centers of stellar photospheric absorption lines 
can also be used, but while these features are located at 
accessible wavelengths, they are broad and shallow and 
therefore hard to identify at the S/N of most high redshift 
galaxy spectra. If none of these features are available, the 
systemic redshift can be estimated from the wavelengths 
of the interstellar absorption lines. 

In this study we use the first and third of the above 
methods. For 51 of the 96 galaxies in the sample, strong 
[O II] AA3727, 3729 emission is included at the red end of 
the spectral region covered, allowing a precise determi- 
nation of the systemic redshift from a double Gaussian 
fit to the [O II] doublet, which is usually resolved. These 
51 galaxies have a mean redshift (z) — 1.399 ± 0.20. 

For the remaining 45 galaxies in the sample, we turn 
to the interstellar absorption lines. It is well-documented 
that these features are blueshifted with respect to the 
systemic redshift and that, while the degree of blueshift- 
ing varies for individual galaxies, the assumption of an 
average velocity offset between the systemic and inter- 
stellar absorption redshifts allows the construction of 
com posite spectra which are effectively at zero veloc- 
ity (IShaplev et alJ 12003b ISteidel et al.ll20lol ). We follow 
the method of ISteidel et alJ (|2010f ) . who use a sample 
of star- forming galaxies at z ~ 2.2 with both interstellar 
absorption lines and systemic redshifts from Ha emission 
to find that, on average, the interstellar absorption lines 
are offset from the systemic redshift by Av = — 166 ± 130 
km s^ 1 . 

We use the subsample of galaxies with measurements 
of both [O II] emission and Fe II A2374 absorption to 
derive the average offset between nebular emission and 
interstellar absorption for our sample (as discussed in 
Sections [5] and [6j this Fe II transition is expected to be 
least affected by the presence of emission, and therefore 
most representative of the kinematics of the outflow) . As 
described in more detail in Section^ absorption line cen- 
troids are measured from a single Gaussian fit using the 
IDL routine mpfit.pro, which measures the center and 
equivalent width of each feature and their associated un- 



certainties. Because [O II] and this Fe II transition are at 
opposite ends of the accessible wavelength range, there 
are 20 galaxies whose spectra include significant detec- 
tions of both features; the average uncertainty in the Fe II 
velocity measurements of this sample is 41 km s _1 . The 
relationship between the [O II]-defined systemic redshift 
and the redshift of the interstellar absorption is given by 

Zsys — ^-abs ( 1 ~l~ ^-abs ) 3 ( 1 ) 

where Ai> a b s = —85 km s _1 is the average velocity offset 
(we compute an unweighted mean to avoid the possibil- 
ity of biasing the result toward brighter galaxies, which 
may also be more massive). The standard deviation is 
87 km s _1 , reflecting the significant scatter among indi- 
vidual measurements, while the error in the mean is 10 
km s _1 . 

For the higher S/N spectra with systemic redshifts 
from [O II] emission, we also measure velocity offsets 
through a joint fit to the F e II A2344, 2374 an d 2587 
lines (for more details, see iMartin et al.l I2012D . The 
mean absorption line velocity measured in this way is 
Av = —76 ± 6 km s _1 with a standard deviation of 80 
km s _1 , in good agreement with the value of —85 km s _1 
measured from Fe II A2374 alone; the joint and Fe II 
A2374 velocities agree within 2a in all cases, and within 
la in half of the measured sample. 

The magnitude of the velocity offset Aw a b s is a fac- 
tor of ~ 2 smaller than the value of Av = —166 ± 130 
km s -1 (mean and standard deviation) determined by 
ISteidel eta l. (2010). We plot Aw a bs against redshift for 
th e current sample an d for the higher redshift sample 
of ISteidel et al l (120 10) in Figure [2] (see also Figure 2 
of ISteidel et al.ll2010() . clearly showing the smaller offset 
from systemic velocity in the lower redshift samp l e. The 
horizontal dotted line divides the [Stcidcl ct "aLl (|2010l ) 
sample into two redshift bins containing an equal num- 
ber of galaxies, and the large light blue points show the 
mean velocity offset and redshift of each bin; although 
the scatter is considerable, the average offset from sys- 
temic increases with increasing redshift, from At; = —110 
km s" 1 at z ~ 2.1 to Av = -203 km s _1 at z ~ 2.4. This 
increase is apparently driven by a decreasing fraction of 
objects with velocities near zero, as well as a higher frac- 
tion with Av < —300 km s _1 . 

Our lower redshift sample follows the same trend of de- 
creasing velocity ccntroid with decreasing redshift. A full 
explanation of this trend requires more careful account- 
ing of selection effects, however; our galaxies are on aver - 
age somewhat lower in mass than the lSteidel et al.l (j2010f l 
sample, since the magnitude limit of the survey reaches 
intrinsically f ainter objects at low er redshifts. Velocity 
offsets in the Stei del et al.l (|2010f ) sample are measured 
from the far-UV interstellar absorption lines C II A1334, 
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Fig. 2. — Velocity centroids of interstellar absorption lines in 
galaxies at z ~ 2 (blue circles, far-UV lines fromfStcidcl ct al.120101) 
and z ~ 1.5 (the current Fe II A2374 sample, green triangles with 
error bars) plotted against redshift. The velocity of Lya emission 
is also shown when present for the z ~ 2 sample (open red squares) . 
The horizonal dotted line divides the z ~ 2 sample into two redshift 
bins containing an equal number of galaxies, and the larger, lighter 
colored points show the mean velocity and redshift of the two bins 
of the z ~ 2 sample and the z ~ 1.5 sample. Velocities are relative 
to Ha emission at z ~ 2 and [O II] emission at z ~ 1.5. 

Si IV A1393 and Si II A1526, whereas we use the near-UV 
Fe II lines. The kinematics of these different lines have 
not been directly compared, although C II, Si II and Fe II 
are expected to be similar as all of these t r ansiti ons arise 
in cool outflowing gas, and lShaplev et al.1 (|2003|) find the 
kinematics of Si IV to be similar to those of the low ion- 
ization lines. Finally, a decrease in the blueshift of the 
velocity centroid does not necessarily correspond to a de- 
crease in the outflow velocity, since the centroid may be 
influenced by the interstellar medium of the galaxy or by 
underlying emission near zero velocity. 

For the 45 galaxies in the sample without [O II] emis- 
sion, we use Equation [T] to determine the systemic red- 
shift. Where possible, we use Fe II A2374 to define the 
absorption redshift, and when that feature is unavailable, 
we use Fe II A2587 (the choice of these two lines is moti- 
vated by the fact that they are least likely to be affected 
by emission filling, which may shift the wavelength to 
the blue). In the 14 galaxies where these two features 
were not detected, we use the far-U V absorption redshif t 
from prior Keck/LRIS spectroscopy (jSteidel et al.ll2004D . 
The mean systemic redshift of the sample of 45 galaxies 
determined in this way is then (z) — 1.836 ± 0.16. Given 
the possible trend of Au a b s with redshift, there is some 
concern that the galaxies for which we apply Equation [T] 
lie at somewhat higher redshifts than the galaxies used 
to determine the equation; however, the similar veloc- 
ity offsets of our samp le and the lower redshift half of 
the iSteidel et al.l i|2()l()i sample suggest that this is not 
a large effect. 

It is desirable to verify systemic redshifts obtained 
in this way via measurement of weak stellar or nebu- 
lar features in composite spectra; if systemic redshifts 
have been assigned correctly, these feat ures should lie at 
zero velocity (e.g. Sh aplev et al.l 120031 ) . Unfortunately 
our spectra do not contain features ideal for this pur- 



pose. The strongest nebular emission line, C II] A2326, 
is a blend of several closely separated transitions whose 
ratios depend on the electron density, making it unsuit- 
able for precise redshift determinations. The composite 
spectrum of all 96 galaxies does contain a weak detection 
of the C III A2297 stellar photospheric line (vacuum rest 
wavelength 2297.579 A), used bv lPettini et~all (|2002[ ) to 
measure the systemic redshift of the lensed galaxy MS 
1512-cB58. The low S/N of this line makes it of lim- 
ited use, but its velocity in the composite spectrum is 
Av = 19 ± 78 km s _1 , consistent with zero. We can also 
check the systemic redshifts by comparing the composite 
spectrum of the 51 galaxies with [O II]-based redshifts 
with that of the 45 galaxies with redshifts determined 
using Equation Q] For this test we compare the centroid 
and width of Fe II* A2626, the strongest emission line ap- 
pearing in both spectra; we find that the line centroids 
in the two spectra agree to within 20 km s , while the 
line is ~ 0.8 A (90 km s _1 ) broader in the spectrum con- 
structed using redshifts from Equation [T] Such broaden- 
ing is expected given the 87 km s _1 dispersion in Au a b s - 
These tests validate our use of Equation Q] to determine 
the systemic redshifts of galaxies without [O II] emission. 

4. STELLAR POPULATION MODELING 

For 71 of the 96 galaxies in the sample, we deter- 
mine the stellar mass, age, extinction and star forma- 
tion rate through modeling of the broadband spectral 
energy distribution (SED). We require photometry in 
the K s and/or Spitzer IRAC 3.6 or 4.5 /im bands in 
order to perform the modeling, and most of the re- 
maining 25 galaxies are not covered by this IR imaging; 
we do not expect these 25 galaxies to be significantly 
different from the rest of the sam ple. We model the 
SEP as described bv lShaplev et al.l (120051) andlErb et all 
(|2006bl) ; briefly, we use the iBruzual fe Charlotl (|2003f ) 
stellar popu lation models, a nd as sume constant star for- 
mation, the iCalzetti et all (|2000fl extinction law, and a 
IChabrierl ([20030 initial mass function. Uncertainties in 
the fitted parameters are determined by Monte Carlo 
simulations which perturb the input photometry accord- 
ing to the photometric uncertainties and determine the 
best-fit SED for the perturbed colors; typical fractional 
uncertainties are (J x /(x) — 0.6, 0.4, 0.5, and 0.2 in 
E(B — V), age, SFR, and stellar mass, respectively. Sys- 
tematic uncertainties are discussed further in Section 
16.11 We assume constant star formation histories be- 
cause of the difficulty of obtaining robust constraints on 
the star formation history of galaxies from SED fitting 
alone; fortunately, the stellar mass is relati vely insen- 
sitive t o the assumed sta r formation history (Er b et all 
200613 IReddv et al.ll20Tl . 

The range of stellar population parameters for the 71 
galaxies is shown in Figure [3] We find a mean (median) 
stellar mass (M*) = 1.1 x 10 10 M (6.3 x 10 9 M ), mean 
(median) reddening (E(B — V)) — 0.21 (0.21), mean (me- 
dian) age 550 Myr (360 Myr), and mean (median) SFR 
54 M Q yr _1 (22 M Q yr _1 ). The degeneracies between 
age, extinction and star formation history involved in 
such modeling are well-known, as is the fact that the 
stellar mass is the most well-determined p a rameter (e.g. 
i Papovich et~aLl l200"lt iShaplev et all 120051: IReddv et all 
120121) . We discuss the relationships between the mod- 
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Fig. 3. — Distributions of stellar population parameters for the 
71 galaxies with SED fits are shown by the pale grey histograms. 
We also divide the sample by redshift and show the distributions 
for the upper and lower halves: the median redshift is z me< j = 1.65, 
galaxies with z < 1.65 are shown with blue histograms, and galax- 
ies with z > 1.65 in red. Galaxies with smaller stellar masses and 
lower SFRs tend to be at lower redshifts because the magnitude- 
limited survey reaches intrinsically fainter objects at lower red- 
shifts. Upper left panel: stellar mass; upper right: age; lower left: 
star formation rate; lower right: E(B — V). 

eled quantities further in Section [6j where we construct 
composite spectra based on stellar population properties. 

5. RESULTS FROM INDIVIDUAL SPECTRA 

Representative examples of the individual spectra are 
shown in Figures [4] and [5] Figure [4] shows the two com- 
plexes of Fe II absorption and emission lines, at ~ 2300 
and ~ 2600 A; the top two spectra are among the highest 
S/N in the sample, while the spectrum of Q2343-BX410 
in the bottom panel is more representative of the typical 
data quality. We also show the Mg II AA2796, 2803 pro- 
files for these four galaxies and four additional objects 
in Figure [5] While the Fe II absorption lines show some 
variation in their strengths, we see much greater diver- 
sity in the Mg II profiles, which range from strong emis- 
sion with very little absorption (Q2343-BX256) to pure 
absorption (Q2343-BX410); the pattern is suggestive of 
that seen in the far-UV spectra of high redshift galax- 
ies, in which Lya may appear in emission or absorption 
while the interstellar resonance lines appear in absorp- 
tion. We study the relationships between the strengths 
of these lines and the properties of the galaxies in this 
section and in the section following. 

The features covered by the spectra are listed in Ta- 
ble O Many of these, particularly the weaker emission 
lines, are significantly detected only in the composite 
spectra discussed in the next section. We measure the 
absorption and emission features by fitting a single Gaus- 
sian profile to each line using the IDL routine mpfit.pro; 
given a spectrum and the error spectrum containing the 
pixel-to-pixel uncertainties, this routine measures the 
center and equivalent width of each line and their associ- 
ated uncertainties. Although absorption lines associated 
with gala ctic outflows generally show asymmetric pro - 
files (e.g. iPettini et all l2000t iQuider et all 120091 [2010h . 



the S/N of the individual spectra in our sample is such 
that weak blue wings of the absorption lines are generally 
undetected and a Gaussian provides a good fit. We do 
not attempt to fit outflowing and zero-velocity compo- 
nents separately, a s some similar stud ies have done (e.g. 
IWeiner et all [2001 iRubin et al.|[2010h . 

5.1. Rest-frame UV Color and Mg II Emission 

We have seen from the spectra shown in Figure [3] that 
the galaxies in our sample show a wide range in Mg II 
AA2796, 2803 profiles, from strong emission to pure ab- 
sorption. In order to separate the sample into galaxies 
with Mg II emission and galaxies without, we classify as 
Mg II emitters those galaxies with at least two adjacent 
pixels at least 1.5a above the continuum level in either of 
the two Mg II transitions. With this classification, 33 of 
the 96 galaxies in the sample show Mg II emission, while 
63 do not. The classification clearly depends on the S/N 
of the individual spectra as well as the intrinsic strength 
of Mg II emission. 

In Figure[5]we plot the rest-frame UV continuum slope 
(i (fx oc A 73 ) against M\jy, the absolute magnitude at 
1900 A, for galaxies with and without Mg II emission, 
using larger symbols for galaxies with higher S /N in the 
Mg II region of the spectrum. We measure f3 from the 
G — 1Z color, which is uncorrelated with redshift for the 
galaxies in our sample. At the sample's mean redshift of 
(z) = 1.6, the G filter is centered at a rest wavelength 
of 1860 A, and the K filter lies at 2665 A. We calculate 
the absolute magnitude Afuv at ~ 1900 A, using the G 
magnitude for 1.4 < z < 1.8, and composite magnitudes 
muG an d mciz for z < 1.4 and z > 1.8 respectively. 
Composite magnitudes arc determined by averaging the 
flux in the two relevant filters. 

It is clear that galaxies with significant Mg II emission 
tend to have bluer UV slopes. Quantitatively, 17 of the 
33 galaxies with Mg II emission are found in the bluest 
third of the sample, with f3 < —1.7, while the remain- 
ing 16 are equally distributed between the two redder 
thirds. A two-sample K-S test finds that the probability 
that the Mg II emitters and non-emitters are drawn from 
the same (3 distribution is 0.025. This conclusion is some- 
what strengthened when we consider only the 52 galaxies 
with S/N > 3 in the Mg II region; among this subsample, 
the 21 Mg II emitters have ((3) = -1.69 ± 0.35 (where 
0.35 is the standard deviation), while the 31 non-emitters 
have (P) = —1.40 ± 0.27. Histograms of these two sub- 
samples are shown in the inset at the upper right of Fig- 
ure [5] Repeating the K-S test indicates that the prob- 
ability that the two samples are drawn from the same 
distribution is 0.0097. In contrast to the (3 distributions, 
the emitting and non-emitting samples are similarly dis- 
tributed in UV luminosity. 

Mg II emission was also identified in a sa mple of ~ 1400 
z ~ 1.4 galaxies by IWeiner et all ([2009), 

using criteria 

which select an excess of emission above the continuum 
level in a window between the two Mg II absorption lines. 
The excess emission objects selected in this way consti- 
tuted ~ 4 % of the sample , and w ere attributed to narrow- 
line AGN. IWeiner et ail ([2009D found that these excess 
emission objects had bluer U — B colors and (in contrast 
to our results) were more luminous. They also found that 
the excess emission was usually at the systemic redshift 



GALACTIC OUTFLOWS IN ABSORPTION AND EMISSION 



7 



1.5 
1 

0.5 


1.5 
1 

x 0.5 

3 

E o 

<D 
N 

£ 1.5 
o 

* 1 
0.5 




_ , 1 1 1 

1 , 


■ ■ ■ ■ ■ 
■ 


,i 


i i i 


i 


. i i | i i . i i | i i 


1 1 

i 


i 


Li i i i 1 i i 


_ : V 




\\ 

III 












^U,l 1 | 1 1 1 1 
h 


- h- - - - - 






i i r 


- ^ — 


i i +■ i i.i i H i i 


l- 




i i 


2325 
i, i i i 


2350 

i i i i i 


2375 

i ,i i 


2400 

i i i.i 


2575 2600 
i i i i i . i i i i i 


2625 
i i i.i 


1 1 




2325 



2350 



2375 



2400 2575 



2600 



2625 







: 1 1 1 : 


. i i | i i.i i 
=- Q2343-BX344 z= 


| i i.i i | , i , i i _ 
1.529 —_ 














- h- -TV 




i .1 i.i i i.i 


i i "t~ i i 


: ~r - i 


1 — i — i ~i — i — h- rr — i 



2325 2350 2375 2400 2575 



2600 2625 




2325 2350 2375 2400 2575 

Rest Wavelength (A) 



2600 



2625 



Fig. 4. — Representative spectra of four individual galaxies, showing the complexes of Fe II absorption and emission. All four galaxies 
have systemic redshifts from [O II] emission. Features are marked by vertical dotted lines. Left panels, C II] A2326 emission, Fe II A2344 
absorption, Fe II* A2365 emission, Fe II A2374 absorption, Fe II A2383 absorption, and Fe II* A2396 emission. Right panels: Fe II A2587 
and Fe II A2600 absorption, Fe II* A2612, A2626 and A2632 emission. 

5.2. Relative Velocities and Equivalent Widths of 
Absorption and Emission Lines 

As discussed in Section [I] most of the Fe II absorption 
lines are coupled to fine structure transitions, such that a 
photon absorbed at (for example) 2587 A can be emitted 
in the Fe II* A2612 or A2632 emission lines. The excep- 
tion is Fe II A2383; the only allowed transition for a pho- 
ton absorbed in this line is a re turn to the ground state . 
Thus, as discussed in detail by IProchaska et aT] (|201lD . 
one may expect significant rescattered emission to "fill 
in" the absorption in this and other resonant transitions; 
one may also expect the filling to be greater when the 
photon has no other means to return to a lower energy 
level, i.e. for transitions such as Fe II A2383 which are not 
coupled to fine structure emission lines. Evidence for this 
emission filling may be seen in the form of blueshifted ve- 
locity centroids, since the scattered emission is expected 
to occur near systemic velocity, or in a decrease in ab- 



and not always accompanied by absorption. 

Our fracti on of Mg II emitter s is considerably higher 
than that of iWeiner et al.l (|2009() . although the different 
selection criteria make a direct comparison difficult. We 
attribute the emission to resonant scattering in outflow- 
ing gas; as shown in Figure [5J the emission is generally 
redshifted and accompanied by blueshifted absorption, 
the classic P Cygni profile of an outflow. The P Cyngi 
profile is also seen in Lya emission and absorption in 
galaxies at high redshift, and the relationship between 
UV color and Lya emission is similar to the relationship 
between UV color and Mg II emission: galaxies with 
stronge r Lya emission are observed to have bluer UV 
slopes (|Shaplev et al.l 120031: Hones et al.l 120111) . We dis- 
cuss the relationship between Mg II emission and galaxy 
properties further in Section [51 
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Fig. 5. — The Mg II AA 2796, 2803 profiles of eight of the galaxies in the sample, showing the range from strong emission to pure 
absorption. Stellar masses are given when available, and galaxies with strong emission tend to be lower in mass. All eight galaxies have 
systemic redshifts from [O II] emission. 



sorption equivalent width. Either of these effects may in- 
fluence conclusions drawn regarding the kinematics, op- 
tical depth and covering fraction of outflowing gas. The 
remainder of this section is dedicated to assessing the ef- 
fect of this emission filling on measurements of velocity 
centroids and equivalent widths. 

In Figure El we compare the velocities with respect to 
systemic for various transitions, as measured by the cen- 
ter of the Gaussian fit. We plot only galaxies with veloc- 
ity uncertainties a\ v < 100 km s _1 (an absolute uncer- 
tainty is preferable to some threshold in <j^ v / Av in order 
to avoid biasing the sample away from galaxies with Av 
near zero) . Because this measurement relies on a precise 
determination of the systemic redshift for each object, 
we include only galaxies for which we have determined 
z sys from [O II] emission. We plot the velocity centroids 
of Fe II A2383 and Fe II A2587 absorption, Mg II A2796 
absorption, and Fe II* A2626 emission, all against that 
of Fe II A2374. As discussed above, we may expect Fe II 



A2383 to show a larger blueshift than Fe II A2374 or Fe II 
A2587, and this is indeed observed in the upper two pan- 
els of Figure El on average, Fe II A2374 and Fe II A2587 
show identical blueshifts of —81 km s _1 , while the aver- 
age relative shifts of Fe II A2374 and Fe II A2383 are -77 
and —129 km s _1 respectively (the two averages for Fe II 
A2374 are not identical because we include only galaxies 
with tJAu < 100 km s" 1 for both lines in each plot, and 
the two subsamples are slightly different). 

We also compare the relative velocity shifts of Fe II 
A2374 and Mg II A2796, in the lower left panel of Fig- 
ure El Because both Fe II A2383 and Mg II A2796 are 
resonance transitions coupled only to the ground state, 
we might expect emission filling to affect them similarly; 
however, as also suggested by the diversity of Mg II pro- 
files relative to Fe II in Figures |H and [5l the Mg II transi- 
tion shows stronger evidence for the presence of emission. 
The average blueshift of Mg II is —164 km s _1 , compared 
to the -129 km s" 1 of Fe II A2383. In the case of both 
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TABLE 2 

Observed Absorption and Emission Features 





\ a 
A lab 


A a 


fa 
J 


Notes 




(A) 


\° ) 




C II] 


2326 


4.43 x 10 1 




Emission; blend 


[O II] 


2470.97 


2.12 x 10~ 2 




Emission; blended with [O II] A2471 




2471.09 


5.22 x 10~ 2 




Emission; blended with [O II] A2470 


Mgl 


2852.96 


4.91 x 10 s 


1.83 


Absorption 


Mg II 


2796.35 


2.60 x 10 8 


0.6155 


Absorption and Emission 




2803.53 


2.57 x 10 s 


0.3058 


Absorption and Emission 


Fen 


2249.88 




0.00182 


Absorption 




2260.78 




0.00244 


Absorption 




2344.21 


1.73 x 10 s 


0.114 


Absorption 




2365.55 


5.90 x 10 7 


0.0495 


Emission 




2374.46 


4.25 x 10 7 


0.0313 


Absorption 




2382.76 


3.13 x 10 s 


0.320 


Absorption 




2396.36 


2.67 x 10 s 


0.279 


Emission 




2586.65 


8.94 x 10 7 


0.0691 


Absorption 




2600.17 


2.35 x 10 s 


0.239 


Absorption 




2612.11 


1.23 x 10 s 


0.122 


Emission 




2626.45 


3.41 x 10 7 


0.0455 


Emission 




2632.11 


6.21 x 10 7 


0.087 


Emission 



a Vacu um wavelengths, Einstein A-values, 
(2003) when available and from the NIST 
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Fig. 6. — Rest-frame UV slope vs. absolute UV magnitude 
Mrjv for galaxies with Mg II emission (blue circles) and without 
(red squares). The symbol size indicates the signal-to- noise ratio 
of the spectrum in the Mg II region: galaxies plotted with large 
symbols have S/N > 4.5 per pixel, medium-sized symbols indicate 
galaxies with 2.5 < S/N < 4.5, and small symbols show galaxies 
with S/N < 2.5. The large, lighter colored symbols show the av- 
erage of the two samples, with error bars showing the standard 
deviations. The inset at upper right shows the distributions 
of Mg II emitters (blue) and non-emitters (red) for spectra with 
S/N > 3. 



lines, the velocity comparison shows that only a fraction 
of the sample shows evidence of emission underlying the 
absorption; about ~ 1/3 of the Fe II A2383 sample is sig- 
nificantly offset from the Fe II A2374 absorption, while 
the fraction of the sample showing an offset in Mg II 
A2796 is higher at - 2/3. 

Finally, we compare the velocity centroids of Fe II 
A2374 with those of the strongest Fe II fine structure 
emission line, Fe II* A2626. We discuss the origin of the 
fine structure emission lines in more detail in the follow- 



and oscillator strengths from Morton 
Atomic Spectra Database otherwise. 

ing sections; for now, we note that (Rubin et all (|2011l ) 
have proposed that Fe II* emission arises from photon 
scattering in the outflowing gas and should be observed 
at or near systemic velocity, assuming that the outflow is 
approximately isotropic and has little extinction. With 
a mean velocity centroid of Ai>2626 = — 20 km s _1 , our 
observations are in general agreement with this scenario. 

We turn next to a comparison of the equivalent widths 
of the Fe II absorption lines, measured by integrating the 
best-fit Gaussian profiles. For an optically thin transi- 
tion, the equivalent width W of an absorption line de- 
pends on the covering fraction and column density N of 
the ion and the oscillator strength of the transition / 
such that W oc NX 2 f. When the transition becomes op- 
tically thick, the equivalent width depends instead on the 
velocity dispersion and covering fraction of the absorb- 
ing gas, and we then expect (in the absence of emission 
filling) different transitions of the same ion to have the 
same equivalent width. It is therefore instructive to com- 
pare the equivalent widths of some of the observed Fe II 
absorption lines. 

In Figure [5] we compare the equivalent widths of Fe II 
A2374 and Fe II A2383. With oscillator strengths of 
/ = 0.0313 and / = 0.320 respectively, these two lines 
should have very different equivalent widths in the opti- 
cally thin case, with WF e IIA2383/WFeIIA2374 = 10-3. As 
the absorbing gas becomes optically thick, we expect 
W Fe iiA2383/W / FciiA2374 ~ 1. An additional effect is the 
possible presence of emission filling: as discussed above, 
the two transitions differ in that Fe II A2374 is cou- 
pled to a fine structure transition, while Fe II A2383 is 
not, and we therefore expect that Fe II A2383 is more 
likely to be affected by the presence of emission underly- 
ing the absorption line. Such underlying emission would 
lead to a decrease in equivalent width, producing, in the 
optically thick case, otherwise unphysical line ratios of 

cIIA2383/W / FcIIA2374 < 1. 

Figure [8] shows that in most cases 
W / FeiiA2383/W / F ciiA2374 is consistent with unity or 
falls between the saturated and unsaturated limits. 
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Fig. 7. — The velocities of Fe II A2383 absorption (upper left), Fe II A2587 absorption (upper right), Mg II A2796 absorption (lower left) 
and Fe II* A2626 emission (lower right), each plotted against the velocity of Fe II A2374 absorption. Only objects with systemic redshifts 
from [O II] emission are included, and the large, light-colored points show the average value for each line. 

ter while BX364 is not. We also note that the stellar 
masses of Q1623-BM95, Q1700-BM154, Q1700-BX1043 
and Q2343-BX364 are 4.1 x 10 9 M Q , 1.9 x 10 9 M Q , 
3.9 x 10 9 M and 4.5 x 10 9 M Q respectively, all lower 
than the average stellar mass of 1.1 x 10 10 M Q . This 
small sample suggests that galaxies with strong Mg II 
emission are more likely to exhibit emission filling in the 
Fe II transitions as well, and that such emission filling 
may be more likely in low mass galaxies. The fraction 
of galaxies significantly affected by this emission is rela- 
tively small, however, particularly for the Fe II lines. 



However, there are four objects (Q1623-BM95, Q1700- 
BM154, Q1700-BX1043 and Q2343-BX364, shown with 
large red circles in Figure |5]) for which the equivalent 
width ratio of the two lines can only be explained by a 
decreased Fe II A2383 equivalent width due to underlying 
emission. Note that this represents only a lower limit 
on the number of objects affected by emission filling; if 
the intrinsic line ratio falls between the unsaturated and 
saturated limits, as is the case for many of the galaxies 
in the sample, the equivalent width comparison cannot 
distinguish between emission filling and an increase in 
optical depth. There are also several additional objects 
that fall in the unshaded region with lower significance, 
so there may be additional galaxies in the sample with 

W / FcIIA2383/W / FoIIA2374 < 1. 

As shown in Figure^ Q1623-BM95 and Q1700-BM154 
also have strong Mg II emission; Q1700-BX1043 and 
Q2343-BX364 (not shown in Figure [SJ have lower S/N 
in the Mg II region, making the presence of emission 
more difficult to assess, but using the criterion de- 
scribed in Section 15.11 BX1043 is also an Mg II emit- 



6. RESULTS FROM COMPOSITE SPECTRA 

In order to examine weak emission and absorption lines 
and study the variation of spectral features with galaxy 
properties more closely, we construct composite spectra 
based on the modeled stellar population properties and 
on the presence or absence of Mg II emission. We also 
create a composite spectrum of all 96 galaxies in the 
sample. 

The stacked spectra represent pixel-by-pixel averages 
of the normalized individual spectra, created with 3u re- 
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Fig. 8. — A comparison of the equivalent widths of Fe II A2374 
and Fe II A2383. The shaded area shows the allowed region, be- 
tween optically thin and thick transitions; the four marked objects 
lie significantly outside this region, and are discussed in the text. 

jection. The mean and standard deviation of the flux 
values at a given wavelength step were calculated, pixels 
with flux values > 3cr were rejected, and the mean flux 
value of the remaining pixels was output as the final com- 
posite spectrum. Because this was done for each pixel, 
the individual spectra contributing to the mean flux at a 
given pixel vary. The error spectrum for each composite 
was computed by adding in quadrature the error spectra 
for the individual galaxies contributing to the mean flux 
at each wavelength step. The final stacked spectra have 
0.25A per pixel in the rest frame. 

The composite spectrum of the full sample of 96 galax- 
ies is shown in Figure |H1 with the features listed in Table 
[5] marked with vertical black lines. In addition to the 
strong Fe II and Mg II absorption lines, we clearly de- 
tect five Fe II* emission lines and the broad blend of 
nebular C II] emission at 2326 A. Figure [S] also shows 
(overplotted in orange) the synthetic, normalized stellar 
spectrum of an integrated stellar population with con- 
stant star formation, smoothed to a resolution of ~ 180 
km s _1 (FWHM) to match the resolution of our observa- 
tions. The spectrum was produced by the st ellar synthe- 
sis code Starburst99 (|Leitherer et al.l fl999f l . and incor- 
porates the library of theoretical UV spec tra of massive 
stars described by iLeitherer et al.1 (|2010| ). The stellar 
spectrum assumes solar metallicity, an age of 100 Myr 
(the UV spectrum stabilizes after ~ 50 Myr und e r con - 
ditions of constant star formation) and a lKroupal ()2001[ ) 
initial mass function. This spectrum shows that the stel- 
lar contribution to the features of interest for the current 
study is negligible; for constant star formation, the stel- 
lar continuum at near-UV wavelengths is dominated by 
stars too hot to have strong Fe II or Mg II absorption. 

6.1. Composite Spectra Based on Stellar Population 
Properties 

As described in Section |H we model the broadband 
SEDs of 71 of the 96 galaxies in the sample to determine 
their stellar masses, ages, star formation rates and red- 
dening. We divide the 71 galaxies into three bins based 



on each of these four parameters and construct compos- 
ite spectra of the galaxies in the highest and lowest bins, 
resulting in 8 composite spectra, each constructed from 
24 individual spectra. The mean stellar population prop- 
erties of the galaxies in each composite are listed in Table 
G2 along with the mean S /N per pixel of each composite 
spectrum. 

There is, of course, considerable overlap between the 
composites because of covariance in the fitted parame- 
ters. We quantify this overlap in Table |4j in order to 
clarify the relationships between the stellar population 
properties considered. This table shows the number of 
objects in common between any two composite spectra. 
The largest overlap is between mass and age; because the 
models assume constant star formation, the stellar mass 
is simply the product of the age and star formation rate. 
As shown in the table, the mass and age composites have 
18/24 objects in common in the massive, old sample and 
16/24 in the low mass and young sample. 

Recent res u lts i ndicate that the use of the 
iCalzetti et al.l (|2000|) extinction law for galaxies with 
young best-fit ages overestimates their star formation 
rates with respect to other indicators; these galaxies may 
have less extinction than is suggested by their red UV 
slopes and be b etter modeled with a stee per, SMC-like 
extinction law (jReddv et al.l 120101 I2012I ). The result 
is that E(B — V) and the SFR may be overestimated 
for the youngest galaxies in the sample (and because 
of the covariance between extinction and age in the 
modeling, the ages of these objects are also uncertain). 
This systematic uncertainty primarily affects galaxies 
with best-fit ages (calculated with models assuming 
the Calzetti extinction law) less than 50 Myr, of which 
there are 16 i n our sample of 71 objects with SED fits. 
iReddy et~a l. (2012) find that restricting the ages to be 
at least 50 Myr (approximately equal to the dynamical 
time) and using the SMC extinction law for such objects 
brings the SED-determined SFRs into agreement with 
SFRs determined from the combined unobscured UV 
+ IR luminosities. The effect of these changes can be 
significant; the typical result is to decrease the best- fit 
E(B — V) by ~ 0.2 mag, and reduce the SFR by a factor 
of ~ 8. Because of these uncertainties, in the following 
discussions we give greatest weight to the composite 
spectra based on stellar mass. 

We have also constructed composite spectra of the 33 
galaxies with and 63 galaxies without Mg II emission, 
determined as described in Section 15.11 Table [3] does 
not report mean stellar population properties for these 
composites, since not all of the galaxies have photome- 
try sufficient for SED modeling; however, the composites 
with and without Mg II emission have mean S/N per 
pixel of 16. The composite spectra based on the stellar 
population properties and on Mg II emission are shown 
in Figure [TUl 

Equivalent widths and velocity centroids are measured 
by direct integration of the flux above or below the con- 
tinuum level, with ler uncertainties determined by Monte 
Carlo simulations in which we perturb each pixel of the 
original composite spectrum by an amount dependent 
on its uncertainty 5000 times to create 5000 artificial 
spectra. The resulting velocity measurements are flux- 
weighted centroids, measured over a typical range of 
~ 900 km s _1 for the absorption lines and ~ 500 km s _1 
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Fig. 9. — The composite spectrum of all 96 galaxies in the sample. The thin red line shows the la error spectrum, and features listed in 
Table [2] are marked with vertical black lines. The overplotted orange line is the stellar spectrum of a solar metallicity, 100 Myr-old galaxy 
with constant star formation as predicted by Starburst99, as discussed in the text. The mean S/N per pixel of the composite spectrum is 
21. 



TABLE 3 

Mean Stellar Population Properties of Composite Spectra 



Composite 


N 


S/N 


Stellar Mass 


E(B-V) 


SFR 


Age 






(per pixel) 


(10 9 M Q ) 




(M yr" 1 ) 


(Myr) 


High Mass 


24 


14 


24.3 


0.20 


28 


1119 


Low Mass 


24 


11 


1.9 


0.25 


54 


133 


High E(B-V) 


24 


12 


12.0 


0.31 


118 


363 


Low E(B-V) 


24 


15 


7.7 


0.14 


18 


589 


High SFR 


24 


15 


9.5 


0.28 


128 


205 


Low SFR 


24 


11 


8.1 


0.17 


10 


768 


High Age 


24 


11 


21.7 


0.18 


16 


1260 


Low Age 


24 


12 


2.7 


0.28 


118 


55 



TABLE 4 

Overlap Between Composite Spectra a 





Low 


High 


Low 


High 


Low 


High 


Low 


High 




Mass 


Mass 


Age 


Age 


E(B-V) 


E(B-V) 


SFR 


SFR 


Low Mass 






18 


1 


8 


10 


9 


9 


High Mass 









16 


9 


6 


6 


7 


Low Age 


18 









5 


15 


3 


15 


High Age 


1 


16 






12 


5 


13 


2 


Low E(B-V) 


8 


9 


5 


12 






13 


3 


High E(B-V) 


10 


6 


15 


5 






2 


15 


Low SFR 


9 


6 


3 


13 


13 


2 






High SFR 


9 


7 


15 


2 


3 


15 







a Number of galaxies in common between any two composite spectra, 
composite represents the average of 24 galaxies. 



Each 



for the emission lines. The uncertainties adopted for Wq 
and velocity are the standard deviation of the 5000 mea- 
surements of these values made for each feature. Mea- 
surements of the equivalent widths and velocity centroids 
of the absorption features are given in Table [5J and of the 
emission features in Table [6] 

Although we measured an average velocity shift of —85 
km s _1 relative to systemic for the Fe II A2374 absorption 
line and then applied this shift to galaxies without [O II] 
detections in order to construct the composites, the ve- 
locity centroid of Fell A2374 in the composite spectrum 
of all galaxies is —120 km s _1 . This difference is due to 
the asymmetric absorption line profiles of the higher S /N 
composite, in which an extended blue tail of absorption 



shifts the centroids to larger velocities when they are 
measured with direct integration rather than Gaussian 
fits which assume symmetry. The peak optical depth of 
the line is consistent with the —85 km s _1 measurements 
from individual spectra. 

The most striking feature of the composite spectra 
shown in Figure [10] is the difference in Mg II profile, 
particularly for the high and low mass composites. The 
lowest mass galaxies have strong Mg II emission and 
weak absorption, while the most massive galaxies have 
stronger absorption and no apparent emission. Galax- 
ies with young ages, low dust content, and (marginally) 
lower SFRs also show stronger Mg II emission, although 
the differences are not as striking as in the compos- 
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Fig. 10. — Composite spectra constructed as a function of stellar population properties and the presence or absence of Mg II emission. 
The spectra focus on the Fe II transitions at ~ 2350 A in the left panels, the Fe II transitions at ~ 2600 A in the middle panels, and Mg II 
AA2796, 2803 in the right panels. From top to bottom, we show spectra binned by stellar mass, E(B — V), age, star formation rate, and 
Mg II emission. Red lines show the high mass, dust, age and SFR subsets, and the galaxies without Mg II emission. Blue lines show the 
low mass, dust, age and SFR subsets, and the galaxies with Mg II emission. 
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Fig. 11. — Velocity comparisons for composite spectra. The velocities of Fe II A2383 absorption (upper left), Fe II A2587 absorption 
(upper right), Mg II A2796 absorption (lower left) and Fe II* A2626 emission (lower right) are each plotted against the velocity of Fe II 
A2374 absorption. Compare Figure [7] for the same measurements using individual spectra. 
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ites based on stella r mass. Similar trends were seen by 
iWeiner et al.l (|2009ft . who saw increased Mg II emission 
in composite spectra of lowe r mass galaxies and galaxies 
with lower SFRs (although Weincr ct al. 2009| removed 
objects with excess Mg II emission from their sample, ad- 
ditional emission became apparent when the symmetric, 
zero-velocity absorption component was removed from 
the Mg II profile; we do not attempt such decomposi- 
tion here) . The composite spectra also show more subtle 
differences in the strengths of the other absorption and 
emission features. We discuss these differences below. 

6.2. Relative Velocities and Equivalent Widths of 
Absorption and Emission Lines in Composite 
Spectra 

An inspection of the absorption line measurements 
in Table [5] shows that the most consistent differences 
between the composite spectra arise when galaxies are 
divided based on stellar mass. All of the lines are 
blueshifted to larger velocities in the high mass spec- 
trum relative to the low mass, and most of the lines have 
higher equivalent widths in the high mass galaxies as 
well. This larger blueshift suggests that more massive 
galaxies drive faster outflows, although we caution that 
the velocity centroid of the line is a poor proxy for out- 
flow velocity because it may be influenced by absorption 
from the interstellar medium of the galaxy or by underly- 
ing emission. We discuss the question of outflow velocity 
in more detail in Section f6. 2. 31 

6.2.1. Relative Velocities of Absorption and Emission Lines 

In Figure [TT] we compare the velocity centroids of sev- 
eral of the lines in the various composite spectra. This 
figure should be compared with Figure in which we 
made the same comparisons for the individual spectra 
and found some evidence for the presence of underly- 
ing emission (in ~ 1/3 of the sample) in the form of 
blueshifted Fe II A2383 velocities relative to Fe II A2374 
and Fe II A2587; recall that the latter two lines are cou- 
pled to fine structure transitions and therefore may be 
less affected by emission filling. This signature blueshift 
of emission filling is less apparent in the composites. 
We see generally good agreement between all three of 
the Fe II absorption lines, though the scatter is signifi- 
cant, particularly in the comparison of Fell A2374 and 
Fell A2587. We conclude that although underlying emis- 
sion may affect the velocity centroids of some Fe II lines 
in some galaxies, the effect is not strong enough to be 
robustly detected in composite spectra. This may be 
because the fraction of galaxies so affected is relatively 
small, or because the higher S/N of the composites in 
combination with the method of direct integration makes 
the measurements of the composites more sensitive to ab- 
sorption at a wider range of velocities, thereby decreasing 
the effect of underlying emission near zero. 

Measurements of the individual spectra showed that 
the velocity centroid of Mg II A2796 is more affected by 
underlying emission than any of the Fe II lines, and this is 
also true of the composites. In contrast to the Fe II lines 
studied, Mg II is blueshifted with respect to Fe II A2374 
in all of the composite spectra. This shift is probably due 
to underlying emission. The largest velocity difference is 
observed in the composites of galaxies with and without 



Mg II emission: the average blueshift of Mg II A2796 
in galaxies with no detectable Mg II emission is Av — 
— 189 ± 24 km s _1 , while the average shift of those with 
emission is Av = — 337 ± 17 km s _1 . We observed above 
that low mass galaxies show stronger Mg II emission, and 
that the emission is somewhat stronger in galaxies with 
young ages, low dust content, and lower SFRs (all shown 
in blue in Figure [TTj) as well. Figure [Til suggests that this 
increased emission in low mass galaxies is apparent in the 
velocity of Mg II absorption as well, since the blue points 
tend to be more displaced than the red points from the 
line of equal velocities. 

Finally, we compare the velocity centroid of the Fe II* 
A2626 emission line with that of Fe II A2374 in the com- 
posites (we note that the velocity of the Fe II* A2626 
line is not affected by the nearby Fe II* A2632 line be- 
cause the two lines are well-separated, as can be seen 
from Figure fT0|) . The individual spectra showed that the 
velocity of Fe II* A2626 scattered around zero, with a 
marginal tendency for greater blueshifts; in the composi- 
ties we see consistently blueshifted velocity centroids of 
about —50 km s _1 . As an additional check on this line, 
we measure the velocity offset of Fe II* A2626 in the com- 
posite spectrum of the 51 galaxies with [O II] redshifts 
and find it to be consistent with the other composites 
at —43 km s _1 . The difference between the individual 
spectra and the composites may again be due to greater 
sensitivity to weak blueshifted emission in the compos- 
ites; see the discussion of the fine structure emission lines 
in Section [6. 31 b clow . In any case, as with the individual 
spectra, Fe II* A2626 is on average observed at velocities 
closer to systemic than any of the Fe II absorption lines. 

6.2.2. Equivalent Widths of Fe II Absorption Lines 

Next we compare the equivalent widths of the Fe II 
A2374 and Fe II A2383 absorption lines in the com- 
posites, as we did for the individual spectra in Sec- 
tion 15.21 and Figure [51 In the case of the individual 
spectra, we saw that, with the exception of four low- 
mass galaxies, all of the objects in the sample have 
an equivalent width ratio consistent with the range ex- 
pected based on the variation of optical depth, 1 < 
H ; FeiiA2383/W / FciiA2374 < 10.3. The ratio of equivalent 
width for all of the composites falls within this range as 
well, as shown in Figure [P2l We do see systematic differ- 
ences in the equivalent width ratio, however: low mass, 
young, and unreddened galaxies with low SFRs (blue 
points) have WFCIIA2383/WFCIIA2374 ~ 1, while older, 
dustier, more massive galaxies with higher SFRs (red 
points) have WfcIIA2383/WfcIIA2374 > 1- These system- 
atic differences may represent real differences in optical 
depth between the two samples, although it may per- 
haps be surprising that the lower mass and less reddened 
samples would show higher optical depth than the dustier 
and more massive galaxies. The alternative explanation 
is that underlying emission is more likely to decrease the 
equivalent width of Fe II A2383 relative to Fe II A2374 
in the low mass sample. This second explanation is sup- 
ported by the four low mass galaxies whose individual 
spectra require the presence of such underlying emission 
(see Figure [5]) , and perhaps by the fact that Mg II emis- 
sion is clearly stronger in low-mass galaxies. However, 
such underlying emission must have relatively little ef- 
fect on the velocity centroids of the lines, since system- 
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atic differences between the red and blue points are less 
apparent in the velocity comparison of Fell A2374 and 
Fe II A2383 shown in Figure [TT] than in the equivalent 
width comparison in Figure 1121 




Fig. 12.— The equivalent widths of the Fe II A2374 and Fe II 
A2383 absorption lines in the composite spectra, plotted with sym- 
bols given in the key at lower right. The shaded region brack- 
ets the optically thin and optically thick line ratios (because the 
two lines have such different oscillator strengths, the upper limit 
of M / F c nA2383/II ; FcllA2374 = 10.3 does not fall within the region 
plotted). Compare Figure \8\ for the same measurements using 
individual spectra. 



6.2.3. Maximum Velocities of Absorption Lines 

We next examine the maximum blueshifted velocities 
of the absorption lines in composite spectra, a measure- 
ment we could not make for the individual objects due 
to limited S/N. High resolution, high S/N spectra of 
lensed galaxies show that the outflow-related interstel- 
lar absorption lin es extend to typical velocities of —700 
to - 800 km s" 1 dPettini et alJl200Ct [Quider et aD l200l 



IMLOl: iDessauges-Zavadskv et alj|201ofr and similar max- 
imum vel ocities are seen in c o mposite spectra of ga laxies 
at z ~ 2 (jSteidel et alj|2010l) . IWeiner et all ([2009) mea- 
sured the velocity v w % at which Mg II A2796 absorption 
reaches 90% of the continuum level in galaxies at z = 1.4, 
and found it to be correlated with both mass and star 
formation rate. Although measurement of the maximum 
wind velocity is dependent on S/N, it is less likely to 
be affected by the interstellar medium or by underlying 
emission than the velocity centroid measurements dis- 
cussed above. 

We define the maximum velocity as the velocity cor- 
responding to the wavelength at which the absorption 
profile first meets the normalized continuum (i.e., has a 
value of one) on the blue side of the line. Because this 
measurement is highly dependent on S/N, we use the 
5000 Monte Carlo simulations described above, and take 
u m „ to be the average of the velocity at which the line 
meets the continuum in the 5000 realizations, with la 
error given by the standard deviation. We measure v max 
for the five strong Fe II absorption lines and for Mg II 
A2796 (we do not measure Mg II A2803 because the blue 



wing is affected by the Mg II A2796 line). The average 
value of v max for these six lines in the composite spec- 
trum of all 96 galaxies is u max = —730 km s _1 , in good 
agreement with other measurements of galaxies at z ~ 2 . 
As other authors have discussed fe.g. lSteidel et al1l2010f ). 
such a velocity is likely to be higher than the escape ve- 
locity of the galaxy, and it is therefore likely that at least 
some of the gas escapes the halo. 

Figure [13] shows u max for each of the composite spec- 
tra. We have averaged the measurements of v max for the 
six absorption lines, and plotted the result against the 
average value of the parameter on which the composite 
spectrum is based (e.g. the average stellar mass of the 
spectrum of low mass galaxies) . The uncertainties in this 
measurement are large, but the maximum outflow veloc- 
ity shows little dependence on dust content, age, SFR or 
the presence of Mg II emission. We do, however, see a 
weak dependence on stellar mass; this is the only param- 
eter for which the maximum velocities of all six lines are 
higher in one composite relative to the other. In the high 
mass composite, we find v max — — 701±131 km s _1 , while 
in the low mass composite we measure u ma x = — 473±116 
km s _1 , for a difference of Av max = 228 ± 124 km s _1 . 
Thus, we find that the maximum outflow velocity in mas- 
sive galaxies is ~ 1.5 times faster that that in low mass 
galaxies, with 1.8cr significance. For comparison, we find 
Av max = —45 ± 117 km s _1 for the E(B — V) compos- 
ites, Ai; max = 20 ± 110 km s _1 for the SFR composites, 
Aw max = 87 ± 132 km s _1 for the age composites, and 
Aw max = —92 ± 132 km s _1 for the composites based 
on Mg II emission. We also note that all of the com- 
posites have similar S/N of 11-16 per pixel (see Table 
[3]); while the S/N of the high mass sample is somewhat 
higher than that of the low mass sample, this is also the 
case for several other pairs of composites, indicating that 
it is not likely to be the cause of the velocity difference 
observed between the high and low mass composites. 

Using 1400 s pectra of gala xies at z ~ 1.4 from the 
DEEP2 survey. IWeiner et all (|2009l ) found that v w% oc 
M°' 17 and v w% oc SFR ' 3 . We clearly do not see the 
same dependence on star formation rate, as our mea- 
sured values of u max are nearly identical in the high and 
low SFR composites. Given the mean SFR values for the 
two composites (10 and 128 M Q yr _1 ), we would expect 
the high SFR sample to have a velocity 2.1 times higher 
that that of the low SFR sample, a result inconsistent 
with the observed AvhighSFft/AviowSFR — 1-0 =t 0-3 at 
the level of 3.7<r (though given the systematic uncertain- 
ties on the SFR discussed above, it would be desirable to 
repeat this measurement with an independent estimate 
of the SFR). Similarly, in a sam ple of galaxies with M g II 
and Fe II absorption at z ~ 1. iKornei et al.l (|2012f ) find 
only weak (~ la) evidence for a trend between outflow 
velocity and star formation rate. Some insight into the 
lack of correlation may be found by considering galaxies 
in the local universe. Using observations of Na I absorp- 
tion in a sample of ne arby ga l axies spanning a range of 
four decades in SFR, iMartinl (|2005l ) found a scaling of 
v oc SFR 35 , while iRupke et al.l (|2005f) suggest that the 
trend of increasing velocity with SFR flattens for SFRs 
above ~ 10 M yr _1 . If this is the case, high redshift 
samples probing a limited range in SFR may be unlikely 
to detect a relationship with velocity. We note, how- 
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Fig. 13. — The average maximum absorption velocity f ma x in the composite spectra, measured as described in the text from the average 
of the five strong Fe II lines and Mg II A2796. 



ever, that the lWeiner et al.l (]2009f ) sample does not span 
a broad range in SFR; their low SFR bin contains galax- 
ies with SFR < 14 M yr -1 , while their high SFR bins 
consists of galaxies with SFR > 28 M Q yr -1 . 

Turning to our observed mass dependenc e, we find that 
it is e ntirely consistent with the results of iWeiner et al.l 
(2009), however; for the mean stellar masses of the up- 
per and lower mass bins, a dependence of v meix cx M® A7 
would result in the massive galaxies having a maximum 
velocity 1.5 times larger than the low mass galaxies, ex- 
actly the difference between our two measured values. 
Given the uncertainties in our measurements of u ma x, 
however, a variety of mass dependencies would be con- 
sistent with our observations. 

6.3. Fe II* and Mg II Emission 

The remainder of this paper is devoted to the emission 
features in the composite spectra. As listed in Table [2] 
and shown in Figures [9] and I10[ we detect emission from 
five Fe II* fine structure lines, and see significant emission 
in the Mg II AA 2796, 2803 doublet for about 1/3 of the 
galaxies in the sample. Mg II emission is more likely 
to appear in galaxies with blue UV slopes and in lower 
mass galaxies. Galaxies with Mg II emission also have 
stronger Fe II* emission; this can be seen most clearly 
in the composite spectra of galaxies with and without 
Mg II emission shown in the bottom panels of Figure 
[TOl and via the measured equivalent widths of the Fe II* 
lines in these spectra given in Table |6] Given the anti- 
correlation between Mg II emission and stellar mass, it 
is then unsurprising that low mass and low SFR galaxies 
also tend to have stronger Fe II* emission compared to 
their higher mass, higher SFR counterparts. 

Fe II* emission is generally not seen in local star - 
forming and star burst galaxies (|Leitherer et al.l [201 lh . 



These lines have also not been seen in samples of galaxies 
at redshifts similar to t he current sample, du e to spectral 
coverage in the case of IWeiner et al.l (|2009f) and proba- 
bly du e to S/N limitations in the sample of iRubin et al.l 
(2010). The situation is similar to that of the Si II* fine 
structur e lines in the far-UV, which are seen i n galaxies at 
z ~ 2 -4 dS haplev et al.l 2003t lErb et al.l2010tlJones et al.l 
|2011[) but not i n loca l starbursts (jSchwartz et al. 2006). 

IRubin et "all (|2011l ) detect both Fe II* and Mg II emis- 
sion in a starburst galaxy at z = 0.69, and propose that 
both are generated by photon scattering in outflowing 
gas. They also show that the Mg II emission is signif- 
icantly spatially extended compared to the stellar con- 
tinuum. This model may explain the absence of fine 
structure emission in local galaxies, if the slits used to 
observe such objects cover only the central regions of the 
galaxies and capture little em is sion from an ex t ended 
outflow (|Giavalisco et alJl2011f >. iProchaska et al.l (|2011l ) 
use radiative transfer models to predict the strength of 
Fe II and Mg II absorption and emission for a variety of 
outflow models; we compare our results to these models 
below. 

As a first test of the origin of the Fe II* emission, 
we compare the velocity structure of the strongest line, 
Fe II* A2626, with that of the absorption line from which 
the emission is expected to arise, Fe II A2600, using the 
composite spectrum of all 96 galaxies. The compari- 
son is shown in the left panel of Figure Q3J with Fe II* 
A2626 emission in blue and Fe II A2600 absorption in 
red. The absorption line shows the clear asymmetric 
profile of galactic outflows, with a blueshifted centroid 
and tail extending to large velocities. Fe II* A2626 is a 
nearly perfect mirror of the absorption profile in emis- 
sion, suggesting that the two lines do indeed arise from 
the same gas. The Fe II* A2626 profile also indicates 
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the presence of dust; since the Fe II* lines are optically 
thin (we see no absorption at their wavelengths), Fe II* 
emission should reach us from both sides of the outflow, 
rather than only the near side seen in absorption. The 
fact that we see more emission from the blueshifted side 
than the redshifted side is probably due to extinction, 
as photons from the far side of the g alaxy must pass 
throu gh more dust in order to reach us ()Prochaska et al.l 
120111) . As a test of this scenario, we also compare the 
Fe II* A2626 profiles in the high and low E(B — V) com- 
posite spectra, in the right panel of Figure [14] If the 
asymmetric line profile is due to extinction, we would 
expect greater asymmetry in the composite spectrum of 
the dustier galaxies. The lower S/N of the extinction- 
based composites makes a robust comparison difficult, 
but Figure [14] docs show that the red wing of the line ex- 
tends to somewhat higher velocities in the low E(B — V) 
composite. A comparison of the equivalent widths of 
the blue and red sides of the lines also offers support to 
the model in which the line asymmetry is due to extinc- 
tion; the red flux is ~ 55% of the blue flux in the high 
E(B — V) composite, and ~ 65% of the blue flux in the 
low E(B — V) composite. 

In this section we examine the origins of Fe II* and 
Mg II emission in more detail, through photoionization 
modeling (Section I6.3.1[) , a study of the spatial extent 
of the emission (Section 16.3.21) , and a comparison of the 
relative strengths of the Fe II absorption lines and their 
associated Fe II* emission (Section 16.3. 3|) . 

6.3.1. Photoionization Modeling 

In this subsection we use photoionization models to 
study the origi ns of Fe II* and Mg II emission. While 
previous work (jRubin et al.ll201lUProchaska et al.ll201lD 
and the kinematics of the Fe II* A2626 line suggest that 
Fe II* fine structure emission arises in the outflow, other 
sources are also possible. In order to assess the possi- 
bility that the Fe II* emission originates in H II regions, 
we have constructed a suite of photoionization models 
which we use to examine the strength of Fe II* emis- 
sion as a function of metallicity and ionization param- 



eter. The models were cons tructed with version 08.00 
of Cloudy, last described by iFerland et al.l (|1998t ). We 
employ an input ion izing spectrum from Starburst99 
(|Leitherer et al.l [T999D . assuming continuous star forma- 
tion and a metallicity of Z — 0.008, or ~ 0.6Z©. We 
assume an electron density n e — 100 cm -3 , consistent 
with the density inferred from the average value of the 
density-sensitive [O II]A3929/[0 II]A3926 line ratio. (We 
find ([O II]A3929/[0 II]A3926) = 1.4 ±0.1, implying den- 
sities between the low density limit and n e ~ 140 cm .) 
We consider gas phase metallicities of 0.2, 0.4, 0.6, 0.8, 
and 1.0 Zq. We use Cloudy's default H II region abun- 
dance set, and in addition consider models with the Fe 
abundance reduced by factors of 0.8, 0.6, 0.4, and 0.2 for 
each of the metallicities adopted. Each model is charac- 
terized by an ionization parameter U , which is defined 
as the ratio of ionizing photon density to total hydrogen 
density, 

U =tS~ . (2) 

where Q is the rate of ionizing photons (number s _1 , 
given by the Starburst99 model used to construct the in- 
put ionizing spectrum), r is the distance from the ion- 
izing source to the inner edge of the cloud, and n is the 
hydrogen density. The models assume a spherical H II 
region, but the geometry is effectively plane parallel if 
the thickness of the cloud is less than O.lro; this is the 
case for all of our models except those with the highest 
values of the ionization parameter considered, log U ~ 0. 
The models use Cloudy's default stopping criterion, in 
all cases stopping when the temperature falls below 4000 
K. 

Results of the photoionization modeling are shown in 
Figure [T51 We plot the ratio of each of the Fe II* lines to 
[O II] AA3626, 3729 against the logarithm of the dimcn- 
sionless ionization parameter U, which represents the ra- 
tio of ionizing photons to particles. (It is more standard 
to measure the strength of a line of interest against a re- 
combination line such as H/3, but [O II] is the only strong 
nebular emission line for which we have measurements.) 
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Fig. 15. — Colored points show the predicted ratios of the luminosities of the Fe II* fine structure emission lines to [O II] A A 3726, 3729 
emission, plotted against ionization parameter. The O/Fe ratio is the same for all models shown, corresponding to Cloudy's default H II 
region abundance ratio. The points are color-coded by metallicity according to the key in the upper left panel. Grey horizontal lines show 
the measured line ratios in individual spectra. 



The [O II] /Fe II* ratio shows little metallicity depen- 
dence but decreases strongly with increasing ionization 
parameter, driven largely by the decrease in [O II] (and 
corresponding increase in [O III]) as ionization parame- 
ter increases. Grey horizontal lines show the measured 
ratios for individual spectra in our sample. 

It is clear that, in order to match the observed ratios, a 
very high ionization parameter of log U > —2 is required. 
We do not have ionization parameter measurements for 
the galaxies in our sample, but measurements of lensed 
galaxies at z ~ 2 suggest that typical star-forming 
galaxies at these red shifts have —2.8 < \ogU < —2.3 
(Hain line et al.l 120091 ) . higher than values observed in 
most local galaxies but lower than needed to produced 
our observed [O II] /Fe II* ratios. We also note two ad- 
ditional factors. First, the measured line ratios are not 
corrected for dust extinction. This is unlikely to be a 
large effect, given the relatively narrow wavelength sep- 
aration, but the effect of dust would decrease the Fe II* 



lines relative to [O II], making the observed [O II] /Fe II* 
ratios higher than the intrinsic ratios and further increas- 
ing the ionization parameter needed to match the data. 
Second, all of the models shown assume Cloudy's default 
H II region O/Fe ratio, which may not be appropriate 
given that an underabundance of the Fe-peak elements 
relative to a-capture elements (such as O) has been ob- 
served in lensed galaxies at high redsh ifts ([Pettini et al.l 
120021: iDessauges-Zavadskv et aLll2010t ). We have tested 
the effect of a decreased abundance of Fe relative to O 
in the Cloudy models; the effect is to weaken the Fe II* 
emission, increase the [O II] /Fe II* ratio, and thereby 
increase the ionization parameter needed to match the 
observations. In other words, a decreased abundance of 
Fe relative to O shifts all of the colored points in Figure 
1151 upwards. 

In summary, it is unlikely that that the galaxies in our 
sample would have ionization parameters high enough 
to produce the required [O II] /Fe II* ratios, and it is 
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therefore also unlikely that the observed Fe II* emission 
originates in the H II regions of the galaxies. 

Next we examine the question of Mg II. In contrast to 
Fe II* (and resonant Fe II) emissi on, emission from M g II 
is observed in local H II regions (jKinnev et al.lll993|) . so 
it is important to assess the likely sources of this emission 
in high redshift galaxies. In Figure [16] we plot the pre- 
dicted ratios of [O II] to Mg II A2796 and A2803 emission 
against the ionization parameter, color-coding by metal- 
licity as in Figure 1151 The observed ratios are shown by 
the grey horizontal lines; unlike the [O II] /Fe II* ratios, 
they are fully consistent with the expected ionization pa- 
rameters of log U ~ —3 to —2. The photoionization mod- 
eling therefore indicates that, in contrast to the Fe II* 
emission, much of the Mg II emission may arise in H II 
regions. Because Mg II is a resonance line, however, pho- 
tons produced in H II regions then scatter in outflowing 
gas, producing the observed blueshifted absorption and 
redshifted emission profiles. 

We have suggested above (Section [5J that Mg II emis- 
sion in star-forming galaxies at z ~ 1-2 is similar to Lya 
emission at z ~ 2-3, in that it exhibits a complex pro- 
file of combined emission and absorption; we have also 
seen that Mg II emission is more common in low mass 
and less dusty galaxies. The photoionization modeling 
supports this connection, indicating that Mg II emission, 
like that of Lya, arises from the resonant scattering of 
photons produced in H II regions. Thus the produc- 
tion of Mg II photons in H II regions is likely to ac- 
count for the difference between the profiles of Mg II and 
that of Fe II lines not coupled to fine structure tran- 
sitions such as Fe II A2383. For example, the spectra 
of Q1623-BM95 and Q1700-BM154, shown in Figures H 
and 03 show weak to moderate Fe II A2383 absorption 
with no associated emission, while their Mg II profiles 
show strong redshifted emission accompanied by weak 
blueshifted absorption. Under the assumption that emis- 
sion is produced by the scattering of continuum photons 
produced in the outflow, one would expect to see simi- 
lar emission in the Mg II and Fe II A2383 lines; instead 



we see much stronger emission in the Mg II transitions, 
accounted for by the production of Mg II photons in 
H II regions and their subsequent resonant scattering. 
A comparison of the predicted strengths of Mg II A2796 
and Fe II A2383 emission in the photoionization mod- 
els supports this conclusion; for ionization parameters 
\ogU ~ —3 to —2, emission from Mgll A2796 is ~ 7 
times stronger than Fe II A2383 emission. We therefore 
expect scattered emission from H II regions to be much 
more significant for Mg II than for Fe II A2383. 

More detailed photoionization modeling in the future 
may address the question of whether the [O II] and Mg II 
emission a rise from the same parts of the H II regions; ac- 
cording to iProchaska et all (|2011[ ). Mg II emission from 
H II regions arises primarily from recombinations in the 
outer layers. An additional question which cannot be 
addressed by our current models and data is that of ad- 
ditional sources of [O II] emission. [O II] emissi on from 
diffus ed ionized gas is seen in local galaxies (e.g. iMartinl 
Il997f ). and is likely present in higher redshift galaxies as 
well. However, the Cloudy models do not provide pre- 
dictions for such emission or absolute [O II] luminosities 
with which to compare our measurements, and our ob- 
servations do not have the spatial resolution necessary 
to address this question (though see Section [6. 3. 21 below, 
where we show that the [O II] emission has the same 
spatial distribution as the stellar continuum, indicating 
that most of the emission does not arise in an extended 
outflow). Our observed [O II] /Mg II ratios are consis- 
tent with emission from both ions arising in H II regions, 
but other sources of [O II] emission cannot be ruled out. 
Future, more detailed modeling and observations with 
higher sensitivity and spatial resolution will address this 
question. 

To summarize the results of this subsection, the pho- 
toionization modeling indicates different mechanisms for 
the production of Fe II* and Mg II emission. Although 
both ultimately arise in the outflow, Fe II* emission 
arises from the re-emission of continuum photons ab- 
sorbed in the Fe II resonance transitions, while the res- 



GALACTIC OUTFLOWS IN ABSORPTION AND EMISSION 



21 



onant scattering of photons originally produced in H II 
regions accounts for much of the Mg II emission. 

6.3.2. The Spatial Extent of Fe II* and Mg II Emission 

If the Fe II* emission originates in a large-scale galactic 
outflow, it may have a more extended spatial distribution 
than the stellar continuum; similarly, resonant scattering 
of Mg II photons in outflowing gas may produce spa - 
tially extended Mg II emission (Prochaska e t alj 120111) . 
iRubin et all (|2011[ ) study Fe II* and Mg II emission m a 
galaxy at z = 0.69, finding that Mg II emission is broader 
than the stellar cont inuum and extends to distances of 
> 7 kpc. Likewise, ISteidel et al.l (|2011|) use deep nar- 
rowband images to show that diffuse Lya emission sur- 
rounds galaxies at z ~ 2 to distances of ~ 80 kpc; the 
Lya transition is similar to Mg II in its radiative transfer 
properties. 

In order to assess the spatial extent of Mg II and Fe II* 
emission in our sample we create two-dimensional com- 
posite spectra and compare the spatial extent of the line 
emission relative to the stellar continuum. We also con- 
struct a 2D composite [O II] spectrum for comparison. 
The 2D composite spectra are created by cutting out a 
100 A region around the line of interest in each individ- 
ual spectrum. These cutouts are also centered spatially 
(but without spatial resampling) on the object, using 
the object positions stored in the image header by the 
DEEP 2 reduction pipeline. The cutouts are then resam- 
pled to a common dispersion of 20 km s _1 per pixel using 
linear interpolation, and the appropriate 2D wavelength 
solution is applied to each cutout. Because the orig- 
inal two-dimensional spectra are rectified and because 
we consider only a short range in wavelength, the ob- 
ject traces are not curved and the wavelength does not 
depend on spatial location along the slit. We then mea- 
sure the centroid of the object position in each spectrum 
(since the positions in the header are not exact), and use 
these measurements to construct spatial offsets used to 
combine the individual spectra. Further details of the 
combination depend on the line considered and are de- 
scribed below. 

The Mg II composite includes only the 33 galaxies with 
Mg II emission, determined as described in Section 15.11 
Because the peak velocity of Mg II emission varies signif- 
icantly in individual spectra, we shift each 2D spectrum 
both spectrally and spatially, centering on the peak of 
the emission, and combine the individual spectra using a 
S/N- weighted average. We measure the extent of Fe II* 
emission using the strongest of the Fe II* lines, Fe II* 
A2626. We combine the full sample (95 of the 96 galax- 
ies in the sample have spectral coverage of the Fe II* 
A2626 line), aligning each spectrum spatially and again 
combining using a S/N- weighted average. For the [O II] 
composite, we combine the 51 galaxies with [O II] emis- 
sion, again aligning each spectrum spatially and combin- 
ing using a S/N- weighted average. In all three cases we 
then average the two sides of the composites in order 
to increase S/N. We determine uncertainties by creat- 
ing 500 artificial 2D composite spectra using bootstrap 
resampling: we measure the average and standard devia- 
tion of the line and continuum profiles in the 500 artificial 
spectra, and take the la uncertainties in the spatial pro- 
files to be the standard deviation of their measurements 
in the sample of artificial spectra. The two-dimensional 



composite spectra for the three lines are shown in Figure 

El 

Spatial profiles of the line emission and stellar contin- 
uum are shown in Figure 1181 The top two panels show 
Mg II emission: at left, we compare the Mg II emission 
profile (red) with that of the stellar continuum (black), 
where the continuum has been scaled by a factor of 1.6 to 
have the same peak value as the Mg II emission. The blue 
dotted line shows the point spread function, measured 
from the average spatial profiles of two stars appearing 
on the masks. In the upper right panel we show the dif- 
ference between the line and continuum emission. The 
Mg II emission is slightly stronger at larger distances, 
reaching a peak excess of 40% greater than the scaled 
continuum at a radius of 0.8 arcsec. 

In the middle two panels we show the same compari- 
son for the weaker Fe II* emission. In this case we scale 
the stellar continuum by a factor of 1.15 to match the 
peak Fe II* flux, and find that the Fe II* emission is also 
somewhat more extended than the continuum, though 
with a smaller excess; the difference between line and 
continuum peaks at ~ 0.7 arcsec, with a 10% excess in 
the line flux relative to the scaled continuum. Although 
this is a smaller excess than we see in the Mg II emis- 
sion, we detect it with comparable significance because 
of the larger sample of galaxies used in the Fe II* com- 
posite. Because the Fe II* emission lines are stronger in 
galaxies with Mg II emission (see Section I6.3[) , a natu- 
ral additional test would be to measure the excess Fe II* 
emission in galaxies with Mg II emission; unfortunately, 
however, the sample of Mg II emitters is too small to 
enable sufficient S/N to make this measurement. 

The bottom two panels of Figure [T8l show the compar- 
ison of the spatial profile of [O II] emission relative to 
the continuum. [O II] emission is expected to arise pri- 
marily in H II regions, and because it is not a resonance 
line, its spatial distribution should trace that of the star- 
forming regions of the galaxy, with any differences with 
respect to the continuum emission arising from differing 
spatial distributions of the hot young stars producing 
the H II regions and the older stars responsible for the 
continuum (or, at larger radii, the possible generation of 
[O II] emission in the outflow). Figure IT51 shows that the 
[O II] emission indeed follows the spatial distribution of 
the continuum closely. There is no excess of line emis- 
sion at ~ 1 arcsec, as is observed for Mg II and Fe II* 
emission; instead the [O II] emission is slightly weaker 
than the continuum at ~ 0.5 arcsec, although the signif- 
icance of this result is marginal. We conclude from this 
test of the spatial distribution of [O II] emission that the 
observed excesses in Mg II and Fe II are not an artifact 
of the comparison of stronger line emission with weaker 
continuum. 

We have presented evidence that both Mg II and Fe II* 
(but not [O II]) emission are slightly more spatially ex- 
tended than the underlying stellar continuum, with ex- 
cess relative to the continuum detected at the ~ 1.5-2tr 
level. This supports the conclusion that these emission 
lines arise in outflowing gas, though higher S/N mea- 
surements are clearly needed. The S/N of the stacked 
spectra we use to make this measurement is only just 
high enough to detect a difference in the spatial distri- 
bution of lines and continuum, but not high enough to 
measure line emission to significantly larger radii than 
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Fig. 17. — Two-dimensional composite spectra of Mg II emission (top), Fe II* A2626 emission (middle), and [O II] emission (bottom). 
The red vertical lines show the wavelength range over which the spatial profiles are measured. 



the continuum. Detecti ons of extended Ly q emission to 
distances of ~ 80 kpc ( Stri de! et al.l 120111 ) suggest that 
such emission could be present at much larger distances, 
but the weakness of Mg II and Fe II* emission relative to 
Lya mean that such detections at high redshift are likely 
out of the reach of current technology. 

6.3.3. Strengths of Fe II* Emission and Fe II Absorption 

Lines 

We conclude our discussion of the origins of Fe II* 
emission with a comparison of the relative strengths of 
the Fe II* emission lines and their associated Fe II ab- 
sorption. As shown in Figure [TJ each of the fine struc- 
ture emission lines is coupled to an absorption line which 
provides the source of the emitted photons. The generic 
expectation for the model in which Fe II* emission is 
produced by photons absorbed in Fe II transitions in the 
outflow is that Fe I I* emission will i n crease as Fe II ab- 
sorption increases. iProchaska et al.1 (j 2 1 If ) provide ex- 
tended discussion of this issue, and use radiative transfer 
models to predict the relative strengths of absorption 
and emission lines for a variety of outflow models. For 
a spherically symmetric outflow in the absence of dust, 
conservation of photons requires that the total emission 
and absorption equivalent widths sum to zero; thus, in 
the absence of dust or particular geometric effects, we 
expect the total emission and absorption strengths to be 
approximately the same. Note also that the absorbed 
photons may be re-emitted as emission at their original 



wavelength rather than through one of the fine structure 
transitions, producing a P Cygni-like profile; whether or 
not this is likely to happen depends on the probabili- 
ties of the transitions involved and on the details of the 
galaxy, including the spatial and velocity distribution of 
outflowing gas, the presence of an interstellar medium at 
zero velocity, and the dust content. 

We show the expected relationship between absorption 
and emission in Figure 1191 which shows the predicted 
equivalent widths of the Fe II A2600 a nd Fe II* A2626 
lines i n all of the models considered by IProchask a et all 
(2011). The points represent a wide variety of outflow 
scenarios rather than a range of galaxies within the same 
model; thus they should be taken simply as a generic pre- 
diction of how increasing one line affects the other. The 
left panel compares the strengths of Fe II A2600 absorp- 
tion and Fe II* A2626 emission, and shows that, although 
there is considerable scatter, Fe II* A2626 emission gen- 
erally increases as Fe II A2600 absorption increases. The 
right panel shows that much of the scatter is due to the 
fraction of emission which escapes at 2600 A rather than 
at 2626 A; except in cases with significant extinction or 
particular outflow geometry, the sum of the total emis- 
sion and absorption is approximately zero. The outliers 
(points falling far from the solid line in the right panel 
of Figure 1 1 9|) represent either non-spherically symmetric 
models or models with significant dust content. These 
models show that either dust or geometry alone can cause 
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Fig. 18.— The spatial profiles of Mg II A2796 (top), Fe II* A2626 (middle), and [O II] AA3727, 3729 emission (bottom) compared to 
the continuum. In the left panels the emission lines are shown in red and the continuum in black. The right panels show the continuum 
subtracted from the emission. Note the differing y-axis scales. 
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Fig. 19.— Equivalent wi dths of Fe II A2600 and Fe II* A2626, as 
predicted by the models of Prochaska ct al. (2011). In both panels 
the x axis shows the equivalent width of Fe II A2600 absorption. 
In the left panel, the y axis shows Fe II* A2626 emission, and in 
the right panel the y axis shows the sum of the equivalent widths 
of Fe II* A2626 emission and Fe II A2600 emission. The solid line 
in the right panel shows equal emission and absorption. Most of 
the scatter in the relationship between Fe II A2600 absorption and 
Fe II* A2626 emission is due to the presence of Fe II A2600 emission. 

large deviations from the simple conservation of photons. 

In Figure [50] we plot the equivalent width of each of 
the Fe II* emission lines against its associated absorp- 
tion line, for each of the composite spectra (note that 
because of the varying overlap of the samples used for 
the various composites, as discussed in Section 16.11 the 
points on these plots are not all independent). We fo- 
cus the following discussion primarily on the Fe II A2600 
absorption and Fe II* A2626 emission lines, since Fe II* 
A2626 is the strongest of the Fe II* lines, but (with some 
exceptions noted below) the paired lines follow similar 
trends. 

As already noted, Fe II* emission is stronger in less 
massive and less dusty galaxies and in galaxies with Mg II 
emission. It is also immediately apparent that we see a 
trend of decreasing Fe II* A2626 emission with increasing 
Fe II A2600 absorption, the opposite of what we expect 
based on the model described above. We have not in- 
cluded Fe II A2600 emission in this comparison, since it 
is not significantly (> 3<t) detected in any of the com- 
posites; however, examination of Figure [TU] shows some 
evidence for redshifted P Cygni Fe II A2600 emission in 
the low mass and low E(B — V) composites. These are 
the only spectra in which the emission is detected at 
> 2cr significance; we find W^ 00 = -0.21 ±0.08 A in the 
low mass composite, and W^ 00 = —0.16 ±0.08 A in the 
low E(B — V) composite. These two composites already 
show relatively weak Fe II A2600 absorption and strong 
Fe II* A2626 emission, so including emission at 2600 
A simply increases the imbalance between total emis- 
sion and absorption in these samples. More generally, 
Fe II A2600 emission decreases as the optical depth o f 
the transition increases, and the iProchaska et al.l ((201 ll ) 
models show that little to no emission is seen when an 
ISM component is present in the galaxy. Our obser- 
vations therefore suggest higher optical depth in higher 
mass and dustier galaxies, possibly due to a m ore signif- 
icant interstellar medium. iWeiner et al. (2009) also find 
that symmetric Mg II absorption centered at zero ve- 
loci ty is significantly str onger in more massive galaxies, 
and (Steidel et al.l (f2010h find more absorption at v ~ 
km s _1 in galaxies with higher baryonic (gas + stellar) 
and dynamical mass at z ~ 2. A more significant ISM in 
more massive galaxies does not account for the observed 
decrease in Fe II* emission with increasing Fe II absorp- 



tion, however, since the presence of an ISM increases 
both the absorption and total emission, while shifting 
the emission from the resonance lines to the fine struc- 
ture transitions. 

The fact that the total absorption is stronger than the 
tot al emission in all s pectra may be attributable to dust; 
the IProchaska et al.l (|201lD models indicate that the ef- 
fect of dust is to suppress the emission with relatively 
little effect on the absorption profile. As noted above, 
the Fe II* A2626 emission profile in the composite spec- 
trum of the full sample shows evidence for extinction 
in the lack of redshifted emission relative to blueshiftcd 
emission. Furthermore, extinction may at least in part 
explain the observed trend in JVabs vs. W cm , if dust 
attenuates a larger fraction of the emission in massive 
galaxies. This does not at first appear to be a satisfac- 
tory explanation, since the average value of E(B — V) is 
slightly higher in the low mass sample than in the high 
mas s sam ple. On the other hand, as discussed in Sec- 
tion 16.11 uncertainties in the extinction law for young 
galaxies suggest that E(B — V) for this sample may be 
overestimated, and independent measurements of extinc- 
tion from Spitzer 24/im imaging indicate that, on aver- 
age, dust attenuati on is larger in more massive galaxies 
(|Reddv et al.l [2010) . Even if this is the case, however, 
the m agnitude of the e f fect d oes not appear to be suffi- 
cient. Procha ska et al.l (|2011[) model the effects of dust 
as a function of the optical depth Td us t , finding that flux 
is reduced by a factor of ~ (1 + Tdust) -1 ; in order to 
suppress the emission to the observed degree in the high 
mass and high dust samples, Td us t > 3 is required, while 
our values of E(B — V) (which do not suffer from the 
systematic uncertainties affecting the lower mass sam- 
p le) indicate Td us t ~ 1.5-2 for these samples. 

IProchaska et al.l f|201 lh also consider a dust±ISM 
model. As discussed above, the effect of an ISM com- 
ponent is to increase both the absorption and total 
emission, while shifting the emission from the resonance 
lines to the fine structure transitions; if the ISM is also 
dusty, the emission is further reduced, particularly in 
the resonance lines. However, the dust±ISM model with 
Tdust = 1 still exhibits increased emission relative to 
the fiducial, dust-free model, indicating that the increase 
in both absorption and emission outweighs the effect of 
dust; this is counter to our observations of weaker emis- 
sion in more massive galaxies with more significant ISM. 
It may be that such a model with somewhat higher dust 
content would sufficiently suppress the emission without 
obscuring it entirely; more detailed modeling is needed, 
and some degree of fine tuning of the dust content would 
seem to be required to match the observations. For the 
moment we conclude that dust is an important factor but 
may not be the only explanation for the observed trend. 

The remaining factor is the geometry of the absorb- 
ing and emitting regions. Anisotropic outflows viewed 
from certain directions may produce observed absorp- 
tion to emission ratios that deviate from the expected 
conservation of photons; in particular, a galaxy with a 
bi-conical outflow observed edge-on will show no absorp- 
tion but strong emission. Thus observations of individual 
galaxies at a variety of inclinations could produce a trend 
similar to the observed W a bs vs. W em relationship. This 
is unlikely to be the explanation for the current sample, 
however, since our data points are for composite spec- 
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Fig. 20. — The equivalent widths of Fe II* emission lines (y axis in each panel) compared with the equivalent widths of the Fe II absorption 
lines from which the emission is assumed to arise (x axis in each panel). The points represent measurements of composite spectra, as shown 
in the legend at lower right. 

Si II* emission but weaker Si II absorption than galaxies 
at z ~ 3, the same pattern we observe in the Fe II lines 
in low and high mass (and age, SFR, and dust) galax- 
ies at z ~ 1.5. They interpret this pattern as evidence 
that Si II photons are absorbed and re-emitted at larger 
distances at z ~ 3 than at z ~ 4. 

If slit losses are in part responsible for the weakness of 
the fine structure emission lines, significant absorption 
in the resonance lines must occur at large distances from 
the galaxies. Recent result s indic ate that this is likely to 
be the case. ISteidel et all d20Toh use stacked spectra of 
paired foreground and background galaxies to measure 
the strength of low ionization absorption as a function of 
impact parameter, finding W a b s — 0.4 A at a mean dis- 
tance of 31 kpc in both the Si II A1260 and A1527 lines. 
The strength of Si II A1260 remains the same at 63 kpc, 
and significant abso rption in Lya is seen at > 100 kpc. 
ISteidel et alJ (|2011f ) also see low surface brightness Lya 
emission extending to ~ 80 kpc. In a study of Mg II 
absorption systems and [O II] emission at 0.4 < z < 1.3, 
iMenard et al.l (|201lD conclude that outflows are the pri- 
mary source of Mg II absorption systems, and that out- 
flowing gas can reach radii of ~ 50 kpc. These results in- 
dicate that significant absorption and emission are likely 
to take place beyond the radius subtended by the slit; 
the 1" and 1.2" slits used in our study subtend distances 
of 8-10 kpc for g alaxies at 1 < z < 2. 

The models of iProchaska et ahl (|201lD also assess the 
effect of slit loss. The fraction of absorption and emission 
at large radius depends on the wind model considered; 
for most models, they find that ~ 1" slits are sufficient to 
capture most of the emission in high redshift galaxies. A 
notable exception is th e rad iation-pressure-driven wind 
model of iMurrav et al.l (|2005l ), for which a slit subtend- 
ing ~20 kpc is required to admit most of the scattered 



tra of galaxies selected without regard to inclination or 
morphologyQ 

The last issue we consider is the size of the Fe II*- 
cmitting region relative to the area encompassed by the 
spectroscopic slit. If fine structure emission occurs in 
outflowing gas at large distances from galaxies, the emit- 
ting region may be larger than the region subtended by 
the slit; the broader the extent of the emission, the higher 
the fraction that may fall outside the slit. In this con- 
text, previous observations of low ionization absorption 
and fine structure emission in far-UV spectra at high 
redshift provide useful insights; the resonant Si II ab- 
sorption lines are coupled to non-resonant Si II* emis- 
sion i n the same fashion as the Fe II lines studied in this 
work. Sha pley et alJ (|2003l ) detected Si II* fine structure 
emission accompanying Si II absorption in the composite 
spectrum of Lyman break galaxies at z ~ 3, finding that 
the emission lines were too strong to be produced in the 
H II regions of the galaxies but much weaker than the 
Si II absorption lines. Although they found no fully sat- 
isfactory explanation for the Si II* lines (outflow mod- 
els were deemed inadequate due to the narrowness of 
the Si II* lines relative to Lya), they suggested that the 
weakness of the Si II* lines could be due to the spectro- 
scopic slit subtending only a small port i on of the emit- 
ting region. More recently, Uones et al.l ([20111 ) compare 
the composite spectra of Lyman break galaxies at z ~ 4 
and z ~ 3, and find that galaxies at z ~ 4 have stronger 

6 Among local galaxies, selecting by dust content may introduce 
an inclination bias, since dustier galaxies are more likely to be edge- 
on. However, among the parent sample from which the galaxies 
observed here are drawn E(B — V) does not increase with apparent 
inclination (D.R. Law, private communication), and in any case 
such a bias would result in stronger emission and weaker absorption 
for the dustier subsample, the opposite of what is observed. 
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Fig. 21. — Two-dimensional composite spectra of Fe II* A2626 emission in the high (top panel) and low (bottom panel) mass subsamplcs. 
Spatial profiles are measured over the wavelength region indicated by the red vertical lines. 

This difference between the emission line profiles in 
the high and low mass galaxies can be further clarified 
by comparisons of the Fe II* line emission and the con- 
tinuum for the high and low mass samples. The relevant 
spatial profiles are shown in the lower two panels of Fig- 
ure!^ A comparison of the line and continuum emission 
in the high mass sample shows that the emission line flux 
is ~ 30% higher at a radius of 0.8-1.1 arcsec; this dif- 
ference is not statistically significant, however, given the 
large uncertainties associated with the small high and 
low mass samples. A similar comparison of the line and 
continuum emission in the low mass sample shows that 
the two profiles are the same to within 5%. Finally, we 
compare the continuum profiles of the high and low mass 
samples. Massive galax ies are observed to be larger than 
low mass galaxies (e.g. lLaw et al~ll2012h . and a compar- 
ison of the continuum profiles to the PSF in the lower 
two panels of Figure [55] suggests that this is also the case 
in our sample. More quantitatively, at a radius of 0.8- 
1.1 arcsec the continuum emission from massive galaxies 
is ~ 40% higher than the continuum emission from low 
mass galaxies, but as with the line and continuum com- 
parison above, this difference is not statistically signifi- 
cant. 

Thus the more extended line emission observed for high 
mass galaxies shown in the upper left panel of Figure [25] 
appears to be due to a combination of extended line emis- 
sion in massive galaxies and the larger sizes of massive 
galaxies; neither of these effects are statistically signifi- 
cant alone, but in combination they produce the moder- 
ately more significant difference between the Fe II* pro- 
files in high and low mass galaxies. In summary, the ob- 
servation of more extended emission from the outflow in 
more massive galaxies supports the scenario in which slit 
losses are more significant for more massive galaxies. A 
definitive test requires the detection of emission at much 
larger radii, but the weakness of the fine structure lines 
is likely to preclude this test for some time. 



photons. 

Given the detection of significant resonant absorption 
at large radius in similar galaxies at z ~ 2, we suggest 
that, along with dust, increasing slit losses in massive 
galaxies may explain the observed trend of decreasing 
Fe II* A2626 emission with increasing Fe II A2600 ab- 
sorption. In this scenario, the smaller ratio of emission 
to absorption strength observed in high mass (age, dust, 
SFR) galaxies is an indication that Fe II photons are 
absorbed and re-emitted at greater distances from the 
galaxy than in lower mass and less dusty objects. This 
is consistent with the observation in Section 16.2.31 that 
massive galaxies have higher wind velocities. We also 
find support for this scenario in obse rvations of a sin- 
gle low mass galaxy; lErb et ail (|2010| ) studied Si II ab- 
sorption and Si II* emission in a compact, unreddened 
galaxy at z = 2.3, finding that the emission and absorp- 
tion strengths were roughly equal. In this galaxy, very 
few of the Si II* photons are lost to dust or fall outside 
the slit. 

A prediction of this slit-loss model is that Fe II* A2626 
emission should be more extended in high mass galaxies 
than in low mass galaxies. We test this by construct- 
ing two-dimensional composite spectra of the high and 
low mass samples as described in Section f6.3.21 again de- 
termining uncertainties from the standard deviation of 
the spatial profiles of 500 artificial 2D composite spectra 
created using bootstrap resampling. The 2D spectra are 
shown in Figure [5T1 and the spatial profiles are shown in 
Figure [55] In the top two panels of Figure [55] we show 
the Fe II* A2626 spatial profiles at left and the difference 
between the profiles in the high and low mass samples at 
right. Although the uncertainties are significant, we find 
that at a radius of 1" the flux of the high mass sample 
is nearly 5 times that of the low mass sample; the dif- 
ference in flux at 1" is 0.69 ± 0.36, for a significance of 
1.9a. We also observe excess emission to ~ 1.5", with 
lower significance. 
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Fig. 22. — Top left: Spatial profiles of the Fc II* A2626 emission line in the two-dimensional composite spectra of the high mass (red) and 
low mass (blue) galaxy samples. Top right: The excess flux in the high mass line emission relative to the low mass. Bottom left: Spatial 
profiles of the Fe II* A2626 emission line (in red) and the continuum (in black) in the high mass subsample. Bottom right: Spatial profiles 
of the Fe II* A2626 emission line (in blue) and the continuum (in black) in the low mass subsample. 



A second important test of the model is a comparison 
of the strength of absorption as a function of impact 
parameter for both high and low mass galaxies, using 
pairs of foreground and background galaxies or galaxies 
with background quasars. The scenario presented here 
predicts stronger low ionization metal line absorption at 
larger radii for more massive galaxies. Such results have 
already been found for galaxies asso ciated with Mg I I 
absorption systems at lower redshifts ()Chen et al.l ;2010), 
and the coming years will see such studies extended to 
the redshifts probed by this sample and higher. 

We have thus far focused discussion on the Fe II A2600 
absorption and Fe II* A2626 emission lines, as Fe II* 
A2626 is the strongest of the Fe II* lines. Although we do 
not detect all of the Fe II* lines in all of the composites, 
the other lines in the UV1 multiplet show the same trend 
of stronger emission and weaker absorption in low mass 
galaxies; we measure only upper limits on Fe II* A2612 
and Fe II* A2632 in the high mass composites. The UV2 
multiplet presents a confusing picture, however. In this 
case we observe weaker absorption and stronger emission 



in the high mass and high dust samples, although galax- 
ies with Mg II emission also show weaker absorption and 
stronger emission. The explanation for these mixed re- 
sults is unclear. In the UV3 multiplet we again observe 
weaker absorption in the low mass sample than in the 
high mass, but the Fe II* A2365 line is too weak to lead 
to useful conclusions regarding the emission. As the UV2 
observations partially conflict with our proposed model, 
further study of these lines is needed. 

We have shown that the average velocity profile of 
Fe II* A2626 emission, photoionization modeling, and the 
extent of the emission relative to the continuum support 
the scenario in which fine structure emission is generated 
by photon scattering in outflowing gas. Greater extinc- 
tion and/or a greater spatial extent of outflowing gas in 
massive galaxies is required to explain the observed de- 
crease in Fe II* emission with increasing Fe II absorption. 

7. SUMMARY AND CONCLUSIONS 

We have presented measurements of rest-frame near- 
UV Fe II and Mg II absorption and emission in a sample 
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of 96 star- forming galaxies at 1 < z < 2. The presence of 
large-scale galactic outflows is indicated by the blueshift 
of interstellar absorption lines relative to the systemic 
velocity. Our primary results are as follows. 

1. For 51 of the 96 galaxies, we obtain systemic red- 
shifts from the centroids of strong [O II] AA 3726, 
3729 emission, which falls at the red end of the 
spectral range covered. These galaxies have an av- 
erage redshift (z) = 1.4. Interstellar Fe II absorp- 
tion lines in these galaxies show an average offset 
from the systemic velocity of —85 km s _1 , less than 
the average offset of —166 km s _1 measured for 
galaxies at z ~ 2. We use this offset of —85 km s _1 
to estimate the systemic redshift of the remainder 
of the galaxies in the sample. We also show that the 
average velocity offset decreases from —200 km s _1 
at z ~ 2.4 to —110 km s' 1 at z ~ 2.1. These higher 
redshift samples have a slightly higher average mass 
than our z ~ 1.4 sample, and have velocity offsets 
measured from far-UV interstellar absorption lines 
rather than the near-UV Fe II lines used in this 
study. More careful accounting of these selection 
effects, and of the effects of the interstellar medium 
on the velocity measurements, is therefore required 
to determine whether the decrease in velocity offset 
with decreasing redshift represents a real trend of 
decreasing outflow velocity. 

2. Individual galaxies in the sample show more va- 
riety in the Mg II doublet than in the Fe II res- 
onance lines. Fe II is always seen in absorption, 
while Mg II ranges from strong emission to pure 
absorption. Emission is more common in galaxies 
with bluer UV slopes, echoing far-UV observations 
of galaxies at z ~ 2-4 in which Lya emission is 
stronger in bluer galaxies. Mg II emission is also 
stronger in low mass galaxies, and fine structure 
Fe II* emission is stronger in galaxies with Mg II 
emission. 

3. We compare the velocities and equivalent widths of 
the Fe II and Mg II transitions in individual spec- 
tra, in order to assess the effect of emission filling 
on measurements of outflow velocity. Underlying 
emission appears to shift the velocity centroid of 
the Fe II A2383 line to the blue in - 1/3 of the 
sample, while the Mg II A2796 line shows such a 
shift in ~ 2/3 of the sample. Studies of individ- 
ual objects suggest that emission filling may be a 
stronger effect in low mass galaxies, but composite 
spectra show little difference in velocity centroids 
between the Fe II lines expected to be most and 
least affected by emission filling. Mg II is much 
more clearly affected by underlying emission. 

4. The average maximum outflow velocity (the max- 
imum extent of the blue wing of the interstellar 
absorption lines in the composite spectrum of all 
of the galaxies in the sample) is v max = —730 
km s _1 , in good agreement with other studies of 
star-forming galaxies at z ~ 1.5 - 2. |w max | is 1.5 
times larger in the high mass subset of galaxies 
than in the low mass subset, consistent with previ- 
ous results indicating w w i n d ~ M° 17 . The velocity 



centroids \v\ are also blucshifted to higher velocities 
in the composite of high mass galaxies than in the 
low mass composite. Unlike some previous studies, 
we see no trend between w max and star formation 
rate. 

5. We present the following evidence in support of the 
model in which Fe II* and Mg II emission is gener- 
ated by photon scattering in outflowing gas: 

• The velocity profile of Fe II* A2626 emission in 
the composite spectrum of the full sample of 
galaxies shows a blueshifted tail very similar 
to that of the Fe II A2600 absorption line, from 
which the emission is presumed to arise. 

• Photoionization modeling indicates that very 
high values of the ionization parameter, 
\ogU > —2, are required in order for the 
Fe II* emission to originate in H II regions. 
The ionization parameters of the galaxies in 
the sample are unlikely to be this high, im- 
plying that Fe II* emission is produced by the 
re-emission of continuum photons absorbed in 
the Fe II resonance transitions in outflowing 
gas. 

• In contrast to Fe II* emission, photoion- 
ization modeling shows that the strength of 
Mg II emission is consistent with its produc- 
tion in H II regions. We conclude that Mg II 
photons from H II regions resonantly scat- 
ter in outflowing gas, producing the observed 
blueshifted absorption and rcdshifted emis- 
sion profiles. 

• We compare the spatial extent of Mg II and 
Fe II* emission relative to the continuum, and 
find that both emission lines are moderately 
more extended. Higher S/N measurements 
are needed to more robustly confirm this trend 
and measure the true extent of the emission. 

6. Using composite spectra, we observe a trend of 
decreasing Fe II* emission with increasing Fe II 
absorption. This trend is the opposite of what 
one would naively expect, if the Fe II* photons 
are originally absorbed in the Fe II transitions. 
We attribute this trend to a combination of in- 
creased extinction in more massive galaxies and 
more extended outflows in massive galaxies, lead- 
ing to greater slit losses for Fe II* photons. Two- 
dimensional composite spectra of high and low 
mass galaxies indicate that emission from the out- 
flow is stronger at a radius of ~ 10 kpc in high 
mass galaxies than in low mass galaxies. 

Our results suggest that Mg II and Fe II* emission 
provide a valuable but so far little-utilized tool for the 
study of galactic outflows. In particular, as telescope 
capabilities develop, studies of the spatial extent of these 
lines may allow important measurements of the physical 
extent of outflowing gas. These lines are accessible from 
the ground over a wide range in redshift, from z ~ 0.2 
to z ~ 2, and complement the far-UV lines studied in 
galaxies at z ~ 2 and above. 



GALACTIC OUTFLOWS IN ABSORPTION AND EMISSION 



29 



We thank the referee for a thorough and constructive 
report which significantly improved the paper, Max Pet- 
tini for useful advice throughout the project, Ralf Ko- 
tulla for assistance with the construction of composite 
spectra, and Jay Gallagher, Claus Leitherer, Norm Mur- 
ray, Alice Shapley and Christy Tremonti for helpful and 
interesting discussions. Michael Cooper provided valu- 
able assistance with the data reduction, and Chuck Stei- 
del, Alice Shapley and Naveen Reddy made possible the 
survey from which the targets in this paper are drawn. 
D.K.E. was supported in part by the Spitzer Fellowship 
Program of the National Aeronautics and Space Admin- 



istration under Award No. NAS7-03001 and the Cali- 
fornia Institute of Technology. A.M.Q. was supported 
by the Marshall Scholarship and the National Science 
Foundation Graduate Research Fellowship. C.L.M. and 
A.L.H. are supported by the National Science Founda- 
tion through grants AST-0808161 and AST-1109288 and 
by the David and Lucilc Packard Foundation. The anal- 
ysis pipeline used to reduce the DEIMOS data was devel- 
oped at UC Berkeley with support from NSF grant AST- 
0071048. We wish to extend thanks to those of Hawaiian 
ancestry on whose sacred mountain we are privileged to 
be guests. 



REFERENCES 



Benson, A. J., Bower, R. G., Frenk, C. S., Lacey, C. G., Baugh, 

C. M., & Cole, S. 2003, ApJ, 599, 38 
Bouchc, N.. Hohcnsee, W., Vargas, R., Kacprzak, G. G., Martin, 

C. L., Cooke, J., & Churchill, C. W. 2012a, MNRAS, 3207 
Bouche, N., Murphy, M. T., Peroux, C, Contini, T., Martin, 

C. L., Forstcr Schreiber, N. M., Genzel, R., Lutz, D., Gillessen, 
S., Tacconi, L., Davies, R., & Eisenhauer, F. 2012b, MNRAS, 
419, 2 

Bouchc, N., Murphy, M. T., Peroux, C, Davies, R., Eisenhauer, 

F., Forster Schreiber, N. M., & Tacconi, L. 2007, ApJ, 669, L5 
Brooks, A. M., Governato, F., Booth, C. M., Willman, B., 

Gardner, J. P., Wadsley, J., Stinson, G, & Quinn, T. 2007, 

ApJ, 655, L17 
Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 
Calzctti, D., Armus, L., Bohlin, R. C, Kinney, A. L., Koornnccf, 

J., & Storchi-Bergmann, T. 2000, ApJ, 533, 682 
Chabrier, G. 2003, PASP, 115, 763 

Chen, H.-W., Wild, V., Tinker, J. L., Gauthier, J.-R., Helsby, 

J. E., Shectman, S. A., & Thompson, I. B. 2010, ApJ, 724, L176 
Coil, A. L., Weincr, B. J., Holz, D. E., Cooper, M. C, Yan, R., & 

Aird, J. 2011, ApJ, 743, 46 
Cooper, M. C, Newman, J. A., Davis, M., Finkbcincr, D. P., & 

Gerkc, B. F. 2012, in Astrophysics Source Code Library, record 

ascl: 1203.003, 3003 
Dalcanton, J. J. 2007, ApJ, 658, 941 

Dave, R., Oppenheimer, B. D., & Finlator, K. 2011, MNRAS, 
415, 11 

Dessaugcs-Zavadsky, M., D'Odorico, S., Schaerer, D., Modigliani, 
A., Tapkcn, C, & Vernet, J. 2010, A&A, 510, A26 

Diamond-Stanic, A. M., Moustakas, J., Tremonti, C. A., Coil, 
A. L., Hickox, R. C, Robaina, A. R., Rudnick, G. H., & Sell, 
P. H. 2012, ArXiv c-prints 

Erb, D. K., Pettini, M., Shapley, A. E., Steidel, C. C, Law, 

D. R., & Reddy, N. A. 2010, ApJ, 719, 1168 

Erb, D. K., Shapley, A. E., Pettini, M., Steidel, C. C, Reddy, 
N. A., & Adelberger, K. L. 2006a, ApJ, 644, 813 

Erb, D. K., Steidel, C. C, Shapley, A. E., Pettini, M., Reddy, 
N. A., & Adelberger, K. L. 2006b, ApJ, 646, 107 

Faber, S. M., Phillips, A. C, Kibrick, R. I., Alcott, B., Allen, 
S. L., Burrous, J., Cantrall, T., Clarke, D., Coil, A. L., Cowley, 
D. J., Davis, M., Deich, W. T. S., Dietsch, K., Gilmore, D. K., 
Harper, C. A., Hilyard, D. F., Lewis, J. P., McVeigh, M., 
Newman, J., Osborne, J., Schiavon, R., Stover, R. J., Tucker, 
D., Wallace, V., Wei, M., Wirth, G., & Wright, C. A. 2003, in 
Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Scries, Vol. 4841, Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Scries, ed. 
M. Iye & A. F. M. Moorwood, 1657-1669 

Faucher-Giguere, C.-A., Keres, D., & Ma, C.-P. 2011, MNRAS, 
417, 2982 

Ferland, G. J., Korista, K. T., Verner, D. A., Ferguson, J. W., 
Kingdon, J. B., & Verner, E. M. 1998, PASP, 110, 761 

Finlator, K. & Dave, R. 2007, ArXiv e-prints, 0704.3100 

Gabor, J. M., Dave, R., Oppenheimer, B. D., & Finlator, K. 
2011, MNRAS, 417, 2676 

Garnett, D. R. 2002, ApJ, 581, 1019 



Giavalisco, M., Vanzella, E., Salimbcni, S., Tripp, T. M., 
Dickinson, M., Cassata, P., Renzini, A., Guo, Y., Ferguson, 
H. C, Nonino, M., Cimatti, A., Kurk, J., Mignoli, M., & Tang, 
Y. 2011, ApJ, 743, 95 

Hainlinc, K. N., Shapley, A. E., Kornci, K. A., Pettini, M., 

Buckley-Geer, E., Allam, S. S., & Tucker, D. L. 2009, ApJ, 701, 
52 

Hartigan, P., Morse, J. A., Tumlinson, J., Raymond, J., & 

Heathcote, S. 1999, ApJ, 512, 901 
Heckman, T. M. 2002, in ASP Conf. Ser. 254: Extragalactic Gas 

at Low Redshift, 292-+ 
Heckman, T. M., Armus, L., & Miley, G. K. 1990, ApJS, 74, 833 
Henry, A. L., Turner, J. L., Beck, S. C, Crosthwaite, L. P., & 

Meier, D. S. 2007, AJ, 133, 757 
Hopkins, P. F., Cox, T. J., Keres, D., & Hernquist, L. 2008, 

ApJS, 175, 390 
Jones, T., Stark, D. P., & Ellis, R. S. 2011, ArXiv e-prints 
Kacprzak, G. G, Churchill, C. W., & Nielsen, N. M. 2012, ArXiv 

c-prints, 1205.0245 
Keres, D., Katz, N., Dave, R., Fardal, M., & Weinberg, D. H. 

2009, MNRAS, 396, 2332 
Kinney, A. L., Bohlin, R. C, Calzctti, D., Panagia, N., & Wyse, 

R. F. G. 1993, ApJS, 86, 5 
Kornci, K. A., Shapley, A. E., Martin, C. L., Coil, A. L., Lotz, 

J. M., Schiminovich, D., Bundy, K., & Noeske, K. G. 2012, 

ArXiv e-prints, 1205.0812 
Kroupa, P. 2001, MNRAS, 322, 231 

Law, D. R., Steidel, C. C, Shapley, A. E., Nagy, S. R., Reddy, 

N. A., & Erb, D. K. 2012, ApJ, 745, 85 
Lehnert, M. D., Heckman, T. M., & Weaver, K. A. 1999, ApJ, 

523, 575 

Leitherer, C, Ortiz Otalvaro, P. A., Bresolin, F., Kudritzki, 

R.-P., Lo Faro, B., Pauldrach, A. W. A., Pettini, M., & Rix, 

S. A. 2010, ApJS, 189, 309 
Leitherer, C, Schaerer, D., Goldadcr, J. D., Delgado, R. M. G., 

Robert, C, Kune, D. F., de Mello, D. F., Devost, D., & 

Heckman, T. M. 1999, ApJS, 123, 3 
Leitherer, C, Tremonti, C. A., Heckman, T. M., & Calzctti, D. 

2011, AJ, 141, 37 
Mannucci, F., Cresci, G, Maiolino, R., Marconi, A., Pastorini, 

G., Pozzetti, L., Gnerucci, A., Risaliti, G., Schneider, R., 

Lehnert, M., & Salvati, M. 2009, MNRAS, 398, 1915 
Martin, C. L. 1997, ApJ, 491, 561 
Martin, C. L. 2005, ApJ, 621, 227 
Martin, C. L. & Bouchc, N. 2009, ApJ, 703, 1394 
Martin, C. L., Shapley, A. E., Coil, A. L., Kornei, K. A., Bundy, 

K., Weiner, B. J., Noeske, K. G., & Schiminovich, D. 2012, 

ArXiv e-prints, 1206.5552 
Menard, B., Wild, V., Nestor, D., Quider, A., Zibetti, S., Rao, S., 

& Turnshek, D. 2011, MNRAS, 417, 801 
Morton, D. C. 2003, ApJS, 149, 205 

Murray, N., Menard, B., & Thompson, T. A. 2011, ApJ, 735, 66 
Murray, N., Quataert, E., & Thompson, T. A. 2005, ApJ, 618, 
569 



30 



ERB, QUIDER, HENRY AND MARTIN 



Newman, J. A., Cooper, M. C, Davis, M., Faber, S. M., Coil, 
A. L., Guhathakurta, P., Koo, D. C, Phillips, A. C, Conroy, 

C, Dutton, A. A., Finkbeiner, D. P., Gerke, B. F., Rosario, 

D. J., Weiner, B. J., Willmer, C. N. A., Yan, Ft., Harker, J. J., 
Kassin, S. A., Konidaris, N. P., Lai, K., Madgwick, D. S., 
Noeske, K. C, Wirth, G. D., Connolly, A. J., Kaiser, N., 
Kirby, E. N., Lemaux, B. C, Lin, L., Lotz, J. M., Luppino, 
G. A., Marinoni, C, Matthews, D. J., Metevier, A., k 
Schiavon, R. P. 2012, larXiv:1203."3"T92l 

Oppenheimer, B. D. & Dave, R. 2006, MNRAS, 373, 1265 
Oppenheimer, B. D., Dave, R., Keres, D., Fardal, M., Katz, N., 

Kollmeier, J. A., k Weinberg, D. H. 2010, MNRAS, 406, 2325 
Papovich, C, Dickinson, M., k Ferguson, H. C. 2001, ApJ, 559, 

620 

Peeples, M. S. k Shankar, F. 2011, MNRAS, 417, 2962 
Pettini, M., Rix, S. A., Steidel, C. C, Adelberger, K. L., Hunt, 

M. P., & Shapley, A. E. 2002, ApJ, 569, 742 
Pettini, M., Steidel, C. C, Adelberger, K. L., Dickinson, M., k 

Giavalisco, M. 2000, ApJ, 528, 96 
Prochaska, J. X., Kasen, D., k Rubin, K. 2011, ApJ, 734, 24 
Quider, A. M., Pettini, M., Shapley, A. E., k Steidel, C. C. 2009, 

MNRAS, 398, 1263 
Quider, A. M., Shapley, A. E., Pettini, M., Steidel, C. C, k 

Stark, D. P. 2010, MNRAS, 402, 1467 
Reddy, N. A., Erb, D. K., Pettini, M., Steidel, C. C, k Shapley, 

A. E. 2010, ApJ, 712, 1070 
Reddy, N. A., Pettini, M., Steidel, C. C, Shapley, A. E., Erb, 

D. K., k Law, D. R. 2012, ApJ, 754, 25 
Reddy, N. A., Steidel, C. C, Erb, D. K., Shapley, A. E., k 

Pettini, M. 2006, ApJ, 653, 1004 
Rix, S. A., Pettini, M., Leitherer, C, Bresolin, F., Kudritzki, R., 

k Steidel, C. C. 2004, ApJ, 615, 98 
Rubin, K. H. R., Prochaska, J. X., Menard, B., Murray, N., 

Kasen, D., Koo, D. C, k Phillips, A. C. 2011, ApJ, 728, 55 
Rubin, K. H. R., Weiner, B. J., Koo, D. C, Martin, C. L., 

Prochaska, J. X., Coil, A. L., k Newman, J. A. 2010, ApJ, 719, 

1503 



Rupke, D. S., Veilleux, S., & Sanders, D. B. 2005, ApJS, 160, 115 
Schwartz, C. M., Martin, C. L., Chandar, R., Leitherer, C, 

Heckman, T. M., k Oey, M. S. 2006. ApJ, 646, 858 
Shapley, A. E., Steidel, C. C, Erb, D. K., Reddy, N. A., 

Adelberger, K. L., Pettini, M., Barmby, P., k Huang, J. 2005, 

ApJ, 626, 698 

Shapley, A. E., Steidel, C. C, Pettini, M., k Adelberger, K. L. 

2003, ApJ, 588, 65 

Soto, K. T., Martin, C. L., Prescott, M. K. M., k Armus, L. 

2012, ArXiv e- prints 
Steidel, C. C, Bogosavljevic, M., Shapley, A. E., Kollmeier, J. A., 

Reddy, N. A., Erb, D. K., k Pettini, M. 2011, ApJ, 736, 160 
Steidel, C. C, Erb, D. K., Shapley, A. E., Pettini, M., Reddy, N., 

Bogosavljevic, M., Rudie, G. C, k Rakic, O. 2010, ApJ, 717, 

289 

Steidel, C. C, Shapley, A. E., Pettini, M., Adelberger, K. L., Erb, 

D. K., Reddy, N. A., k Hunt, M. P. 2004, ApJ, 604, 534 
Strickland, D. K. k Stevens, I. R. 2000, MNRAS, 314, 511 
Tremonti, C. A., Heckman, T. M., Kauffmann, C, Brinchmann, 
J., Chariot, S., White, S. D. M., Seibert, M., Peng, E. W., 
Schlegel, D. J., Uomoto, A., Fukugita, M., k Brinkmann, J. 

2004, ApJ, 613, 898 

Tremonti, C. A., Moustakas, J., k Diamond-Stanic, A. M. 2007, 
ApJ, 663, L77 

van de Voort, F., Schaye, J., Altay, G., k Theuns, T. 2012, 

MNRAS, 421, 2809 
Veilleux, S., Cecil, G., k Bland-Hawthorn, J. 2005, ARA&A, 43, 

769 

Weiner, B. J., Coil, A. L., Prochaska, J. X., Newman, J. A., 
Cooper, M. C, Bundy, K., Conselice, C. J., Dutton, A. A., 
Faber, S. M., Koo, D. C, Lotz, J. M., Rieke, G. H., k Rubin, 
K. H. R. 2009, ApJ, 692, 187 



TABLE 5 

Measured Absorption Features in Composite Spectra 3, 





All 

Galaxies 




Low 

Mass b 




High 

Mass b 


Low 
E(F3-V) C 


High 
E(B-V) C 


Low 
SFR d 


High 
SFR d 




Low 
Age° 


High 

Agc° 


No Mgll 
Emission 


Mgll 
Emission 


N 


96 




24 




24 


24 


24 


24 


24 




24 


24 


63 


33 


WFeIIA2249 


0.23 ± 0.07 




< 0.65 





.31 ±0.08 


0.38 ±0.08 


< 0.28 


< 0.54 


< 0.24 




< 0.43 


< 0.50 


0.40 ±0.09 


< 0.18 


WFeIIA2260 


0.36 ± 0.11 




< 0.55 




< 0.41 


0.47 ±0.13 


< 0.56 


0.57 ±0.18 


< 0.45 




< 0.48 


< 0.67 


0.51 ±0.12 


< 0.21 


WfcIIA2344 


2.56 ± 0.14 


2. 


00 ±0.12 


2. 


69 ±0.16 


2.52 ±0.18 


2.28 ±0.17 


2.27 ±0.20 


2.45 ±0.13 


2 


.03 ±0.14 


2.19 ±0.20 


2.60 ±0.18 


2.06 ± 0.14 


M / FeIIA2374 


1.77 ±0.12 


1 


.84 ±0.14 


1 


.53 ±0.11 


1.87 ±0.11 


1.47 ±0.10 


1.76 ±0.12 


1.47 ±0.08 


1 


.70 ±0.10 


1.70 ±0.18 


1.90 ±0.17 


1.25 ± 0.09 


M / PcIIA2382 


2.26 ± 0.17 


1 


98 ±0.20 


2. 


.00 ±0.12 


1.95 ±0.11 


1.73 ±0.11 


1.79 ±0.12 


1.94 ±0.10 


1 


.73 ±0.13 


2.19 ±0.12 


2.45 ± 0.23 


1.56 ±0.11 


WPCIIA2586 


2.31 ± 0.11 


1 


.99 ±0.16 


2. 


.80 ±0.15 


2.11 ±0.10 


2.76 ±0.16 


2.24 ±0.13 


2.48 ±0.11 


2 


.07 ±0.14 


2.70 ±0.19 


2.32 ± 0.11 


2.14 ±0.13 


WpcIIA2600 


2.55 ± 0.17 


2 


.58 ± 0.24 


2 


.89 ±0.23 


2.34 ±0.10 


3.28 ±0.28 


2.63 ±0.20 


2.82 ±0.20 


2 


.94 ±0.27 


2.75 ± 0.26 


2.64 ±0.17 


2.09 ±0.21 


WMgIIA2796 


1.60 ±0.12 





.78 ±0.25 


3 


.14 ±0.15 


1.57 ±0.18 


1.95 ±0.25 


1.25 ±0.23 


1.70 ±0.17 


1 


,01 ±0.21 


2.44 ± 0.33 


2.10 ±0.17 


1.17 ±0.08 


M / MgIIA2803 


1.53 ± 0.06 





.79 ±0.25 


3 


10 ±0.12 


1.63 ±0.09 


1.82 ±0.22 


1.40 ±0.13 


2.22 ±0.16 


1 


.13 ±0.26 


2.22 ±0.16 


1.95 ± 0.11 


1.17 ±0.06 


WmeIA28B2 


0.53 ± 0.09 




< 0.71 





92 ±0.12 


0.49 ±0.14 


< 0.64 


0.53 ±0.16 


0.65 ±0.19 




< 0.60 


0.66 ± 0.12 


0.49 ± 0.16 


0.49 ± 0.13 


«FeIIA2249 


24 ±56 








-2 ±40 


-21 ±40 














82 ±47 




■yFeIIA2260 


-108 ± 75 










-151 ± 72 




-136 ±93 










-99 ± 59 




"EVaIIA2344 


-71 ± 33 




-31 ± 17 




-97 ± 33 


-73 ± 42 


-52 ±31 


-123 ±41 


-55 ±23 




-55 ± 26 


-105 ±42 


-101 ±38 


-32 ± 31 


■yFeIIA2374 


-120 ±27 




-39 ± 26 




-127 ±23 


-149 ± 25 


-82 ±16 


-128 ±22 


-95 ± 13 




-62 ± 18 


-140 ±45 


-119 ±36 


-93 ± 22 


^FeIIA2382 


-118 ±46 




-37 ± 40 




-114 ±20 


-176 ± 25 


-46 ± 20 


-189 ±25 


-108 ±16 




-62 ± 24 


-103 ± 19 


-112 ±52 


-143 ±26 


fFeIIA2586 


-91 ± 19 




-53 ± 24 




-120 ±27 


-67 ± 12 


-116 ± 24 


-48 ± 18 


-148 ±15 




-106 ±21 


-59 ± 32 


-102 ± 17 


-69 ± 21 


^FeIIA2600 


-189 ±39 




-131 ±43 




-176 ±42 


-124 ± 11 


-173 ±47 


-139 ±29 


-163 ±36 




-171 ±48 


-173 ±45 


-156 ±33 


-171 ±49 


^MgIIA2796 


-281 ±23 




-230 ± 55 




-255 ± 15 


-306 ± 32 


-217 ±43 


-270 ± 33 


-245 ± 26 




-238 ± 47 


-254 ± 43 


-189 ±24 


-337 ± 17 


^MgIIA2803 


-203 ± 7 




-162 ±33 




-193 ± 11 


-188 ± 10 


-231 ± 31 


-204 ±17 


-190 ±20 




-164 ±24 


-190 ± 19 


-149 ± 11 


-260 ± 8 


«MgIA2852 


-130 ±35 








-71 ± 27 


-102 ± 62 




-134 ±46 


-159 ±55 






-61 ±23 


-92 ± 54 


-146 ±54 



a Equivalent widths Wo are given in A in the upper half of the table, and velocity centroids v are given in km s _1 in the lower half of the table. 

b Minimum, maximum, mean and median stellar masses for the low mass composite: 1.3 X 10 s , 3.9 X 10 9 , 1.9 X 10 9 , 1.8 X 10 9 Mq respectively; minimum, maximum, 

mean and median stellar masses for the high mass composite: 1.1 X 10 10 , 9.7 X 10 10 , 2.4 X 10 1C) , 1.5 X 10 10 Mq respectively. 

c Minimum, maximum, mean and median values of E(B — V) for the low E(B — V) composite: 0.07, 0.18, 0.14, 0.15 respectively; minimum, maximum, mean and 
median values of E(B - V) for the high E(B - V) composite: 0.23, 0.46, 0.31, 0.29 respectively. 

d Minimum, maximum, mean and median SFRs for the low SFR composite: 5, 17, 10, 9 Mq yr —1 respectively; minimum, maximum, mean and median SFRs for 
the high SFR composite: 31, 628, 128, 56 Mq yr" 1 respectively. 

c Minimum, maximum, mean and median ages for the low age composite: 7, 180, 55, 23 Myr respectively; minimum, maximum, mean and median ages for the high 
age composite: 571, 4750, 1260, 856 Myr respectively. 



TABLE 6 
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All 


Low 


High 


Low 


High 


Low 


High 


Low 


High 


No Mgll 


Mgll 




Galaxies 


Mass b 


Mass b 


E(B-V) C 


E(B-V) C 


SFR d 


SFR d 


Age e 


Age e 


Emission 


Emission 


N 


96 


24 


24 


24 


24 


24 


24 


24 


24 


63 


33 


WCII1A2326 


-0.74 ±0.07 


-1.18 ±0.29 


> -0.46 


-1.02 ±0.08 


> -0.96 


-0.97 ±0.16 


-0.90 ±0.18 


-0.83 ±0.27 


> -0.76 


-0.68 ±0.13 


-0.76 ±0.12 


WFCIIA2365 


> -0.22 


-0.23 ±0.07 


> -0.28 


-0.26 ±0.05 


> -0.32 


-0.35 ±0.07 


-0.30 ±0.10 


-0.23 ±0.07 


> -0.23 


> -0.32 


-0.31 ±0.07 


WFCIIA2396 


-0.40 ± 0.05 


-0.27 ±0.07 


-0.70 ± 0.09 


-0.45 ±0.06 


-0.61 ±0.12 


-0.78 ±0.08 


-0.57 ±0.10 


-0.35 ±0.08 


-0.47 ±0.10 


-0.32 ±0.07 


-0.61 ±0.08 


Wf b IIA2612 


-0.50 ± 0.05 


-0.57 ±0.10 


> -0.28 


-0.59 ±0.09 


-0.47 ±0.15 


-0.24 ±0.07 


> -0.51 


> -0.51 


> -0.26 


-0.33 ±0.07 


-0.84 ±0.10 


WFCIIA2626 


-0.90 ±0.05 


-0.94 ±0.10 


-0.73 ± 0.08 


-1.01 ±0.08 


-0.69 ±0.10 


-1.17 ±0.09 


-0.89 ±0.10 


-0.76 ±0.10 


-0.75 ±0.08 


-0.78 ±0.07 


-1.18 ±0.09 


WFCIIA2632 


-0.19 ±0.05 


-0.39 ±0.07 


> -0.21 


-0.34 ±0.09 


> -0.15 


-0.40 ±0.12 


-0.24 ±0.08 


-0.22 ±0.07 


> -0.30 


> -0.18 


-0.37 ±0.12 


M / MgIIA2796 


> -0.16 


-0.49 ±0.14 


> -0.06 


-0.44 ±0.07 


> -0.11 


> -0.35 


> -0.09 


-0.41 ±0.13 


> -0.10 


> -0.10 


-1.17 ±0.06 


M / M E IIA2803 


> -0.17 


> -0.38 


> -0.06 


> -0.23 


> -0.13 


-0.25 ±0.07 


> -0.12 


> -0.22 


> -0.12 


> -0.13 


-0.71 ±0.07 


«CII]A2326 


59 ±25 


4 ±87 




95 ±24 




71 ±48 


-3 ±59 


6 ±96 




23 ±49 


45 ±41 


^FeIIA2365 




-23 ± 29 




-84 ± 23 




-63 ± 27 


6 ±48 


2 ±28 






-54 ±35 


^FeIIA2396 


-33 ± 20 


-42 ± 28 


-17 ±33 


5 ±20 


-119 ±41 


31 ± 19 


-121 ± 36 


-64 ± 30 


3 ±38 


-51 ±32 


-17 ±28 


^FeIIA2612 


-17 ±15 


-1±24 




-51 ±25 


-34 ± 61 


29 ±28 








-22 ± 28 


-56 ±27 


«FeIIA2626 


-45 ± 8 


-40 ± 14 


-3 ±18 


-29 ± 14 


-45 ± 20 


-56 ± 10 


-31 ±30 


-52 ± 20 


-40 ± 15 


-47 ± 12 


-44 ± 17 


^FeIIA2632 


-15 ± 28 


-10 ±16 




-34 ± 37 




-21 ± 30 


-20 ± 48 


-18 ±26 






17 ± 55 


«MgIIA2796 




1 ± 19 




-56 ± 12 








1 ±20 






-92 ± 5 


'"MgIIA2803 












19 ±20 










-43 ± 17 



a Equivalent widths Wo are given in A in the upper half of the table, and velocity centroids v are given in km s~ 1 in the lower half of the table. Negative equivalent 
widths indicate emission. 

Minimum, maximum, mean and median stellar masses for the low mass composite: 1.3 X 10 8 , 3.9 X 10 9 , 1.9 X 10 9 , 1.8 X 10 9 M respectively; minimum, maximum, 
mean and median stellar masses for the high mass composite: 1.1 X 10 10 , 9.7 X 10 10 , 2.4 X 10 10 , 1.5 X 10 10 Mq respectively. 

c Minimum, maximum, mean and median values of E(B — V) for the low E{B — V) composite: 0.07, 0.18, 0.14, 0.15 respectively; minimum, maximum, mean and median 
values of E(B - V) for the high E(B - V) composite: 0.23, 0.46, 0.31, 0.29 respectively. 

d Minimum, maximum, mean and median SFRs for the low SFR composite: 5, 17, 10, 9 Mq yr — 1 respectively; minimum, maximum, mean and median SFRs for the high 
SFR composite: 31, 628, 128, 56 M yr" 1 respectively. 

e Minimum, maximum, mean and median ages for the low age composite: 7, 180, 55, 23 Myr respectively; minimum, maximum, mean and median ages for the high age 
composite: 571, 4750, 1260, 856 Myr respectively. 



